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(57) A characteristic value identification method and 
an apparatus therefor which can develop each model to 
be integrated into an entire model, similar to a product 
in which individual parts are combined, are provided. A 
functional model of a part is prepared based on a po- 
tential quantity and a flow quantity representing a 
strength and a quantity of energy applied to the part, a 
steady internal characteristic value of the functional 
model in a steady state is identified, and a transient in- 
ternal characteristic value of the functional model in a 
transient state is identified by using the identified steady 
internal characteristic value. Furthermore, thefunctional 
model having the characteristic value identified by such 
an identification apparatus is incorporated into a virtual 
testing system as a virtual prototype, an internal char- 
acteristic value of the virtual prototype is evaluated by 
providing a driving operation condition and an environ- 
ment condition, actual machine test data obtained by the 
driving operation condition and the environment condi- 
tion are compared with the internal characteristic value, 
and a re-identification is performed depending on the 
comparison result if necessary. 
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in a process 21 , a steady functional model is obtained which reproduces a function, performance, and characteristic 
of a steady state from the functional model of the part. This steady functional model turns into a model including only 
a steady characteristic except a transient characteristic. In a process 22, a steady test of a part which forms an iden- 
tification subject is performed to collect steady test data. In a process 23, the steady internal characteristic value is 

5 identified by using the test data to be provided to the functional model by a process 24. 

[0012] The transient identification process 3 includes the following processes according to claims 6 to 9. Firstly in a 
process 31 , a transient test model of a part which forms an identification subject is prepared, and the test is performed 
to collect transient test data. In a process 32, a transient internal characteristic value of the functional model is identified 
by using the test data. As shown in a process 4, the steady internal characteristic value identified in the process 24 is 

to reflected to the transient identification process 3. 

[0013] Hereinafter, along these processes, the test model and the identification method for identifying the functional 
model will be described taking a motor as an example enabling a brief description. 

[0014] Firstly, the functional model of the process 1 common to each claim will be described by taking a DC (direct 
current) motor as an example. 

15 

<Modeling of the motor> 



(1 ) Basic function of the motor 



20 



25 



[0015] A basic function of a motor is that electric energy is inputted, and the energy is converted into a rotation energy 
to be outputted. As for the energy conversion at this time, a voltage V M and a current l M are applied to the electric 
system of the motor whereby an angular velocity w M and a torque T M are outputted from the rotation system. Also, a 
loss characteristic, an accumulation characteristic, and an additional load are included in the motor as an internal load, 
so that these determine the behavior and the loss characteristic of the motor. 

[0016] A torque coefficient Xt [Nm/A] expressed by the ratio of current and torque, and a velocity coefficient [V/ 
(rad/sec)] expressed by the ratio of angular velocity and induced voltage are known between the electric system and 
the rotation system of the DC motor. An ideal motor model which directly converts an input/output state value with 
these coefficients being made a basic function can be expressed by the following equation: 



30 



35 



T 1 
X T = — = — -ZcfrxKT 8 



x tt ■ — =^-z4>xio 



27i a 
1 P, 



Eq.(l) 



2n a 



[0017] In equation (1), P is the number of poles in armature, "a" is the number of armature parallel circuits, Z is the 
number of all conductors, and <|> is all of the magnetic flux per pole. They are characteristics determined by a motor 

to structure. Also, the magnetic flux <j> is determined by the structure in a separately excited type motor which uses a 
permanent magnet for a field in the same way as above, and is influenced by a state value of the electric system in a 
self-excited motor having a field winding. Eq.(1 ) is a basic one common to a DC rotator applied to a DC generator as well . 
[001 8] The torque coefficient Xt and the velocity coefficient of Eq.(1 ) performs a physical unit system transforma- 
tion from current into torque, or from angular velocity into voltage in the same equation. Therefore, these coefficients 

45 will be hereafter called motor constants expressed by M M . As for the relationship between the motor constants and 
the internal load, the electric system and the rotation system are considered to be the motor internal load with the 
coefficients being bordered. A basic form of the motor model combining the relationship of the basic function of the 
motor constant M M and the internal load is shown in Fig.2. 

[0019] In Fig.2, the input/output state value of the electric system is connected to a power supply, and the input/ 
so output state value of the mechanical system is united with an external machine load. Also, as for the flow of the state 

value, the lower flow represents a flow system in which the current l M of the power supply is converted by M M into a 

flow of the output torque T M which drives the external load, while the upper flow represents a potential system in which 

the angular velocity co M received from the external load is converted into the motor induced voltage by M M to turn into 

a flow of the voltage V M returned to the power supply. 
55 [0020] Thus, the mechanical system and the electric system are united by the strength of energy called a potential 

quantity and the quantity of energy called a flow quantity, thereby enabling energy rules of two dimensions, equilibrium 

of forces and the continuity of velocity, to be satisfied. 
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(2) Functional model of the motor 

[0021] An electric circuit diagram of the motor model is shown in Fig.3. 

[0022] In the electric system of Fig.3, R M is a winding resistance, L M is an inductance, and R c is an insulation re- 
5 sistance of the winding. Also, the rotation system has a moment of inertia J M , a viscous resistance D M , and a friction 

torque T MF . The functional model expressing these relationships is shown in Figs.4(1) and 4(2). It is to be noted that 

such a functional model itself has already been disclosed in the Japanese Patent Application Laid-open No. 9-91334. 

[0023] In Fig.4(1), the right side of the motor constant M M forms the rotation system internal load, and the left side 

thereof forms the electric system internal load. This model has a voltage drop E MB of a brush which supplies the current 
10 i M to the armature, and a friction torque T MP of the rotation system as the additional load. When the influence of both 

is neglected, T MF and E MB can be assumed to be 0. It is to be noted that x M and x L in Fig.4(1) are internal state values. 

Also, the following velocity correction coefficient 8 M is added to the motor constant M M by taking into account the 

variation of the torque coefficient xt and the velocity coefficient . 
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[0024] A government equation of such a motor functional model in Fig.4(1 ) is given by the following equation: 
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[0025] In equation (3), the 1st-2nd rows form state equations, and the 2nd-4th rows form output equations. 
[0026] It is to be noted that Eq.(3) can be derived from the following process shown in Fig.4(1 ); 

• A white circle adds the inputted state values to make the output state value; 

• At a black circle, the outputted state value is equal to the inputted state value; 

• A characteristic expressed by a square makes the product of the input state value and the characteristic the output 
state value; 

• A white circle with a mark x outputs the product of the input state values or input signals; 

• A white circle with a mark - makes the state value negative; 

• A triangle integrates a differential quantity of the inputted state value to make the output state value; 

• A large white circle on which a variable name is written expresses the state value generated in the functional model 
as an additional load. 

(3) Uniting with motor power supply . 

[0027] A battery forming the power supply of the motor is modeled. Fig.5 shows a circuit diagram of the battery, and 
Fig.6 shows its functional model. In the battery of Fig.6, a voltage V B and a current l B are connected to the voltage V M 
and the current l M of the motor shown in Fig.4(1). Also, the characteristic of Fig.6 has E 0 as an induced electromotive 
force of the battery, and R E as an internal resistance of the battery. 
[0028] The government equation of the battery in Fig.5 is given as follows: 
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...Eq.(4) 



[0029] Uniting the motor model and the battery model is given by the following equation by substituting equation (4) 
for Eq.(3). 
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...Eq.(5) 



[0030] The 1st-2nd rows of equation (5) form the state equations, and the 3rd-5th rows form the output equations. 
An abstracted equivalent admittance Y 0 and a current source l 0 in Eq.(5) are expressed by the following equations: 



■o = Y o E o Eq.(7) 

[0031] Thus, a characteristic value identification method according to the present invention uses a functional model 
of a part based on a potential quantity and a flow quantity representing energy applied to the part, and comprises a 
first process for preparing that functional model, a second process for converting the functional model into a steady 
functional model in a steady state to identify a steady internal characteristic value, and a third process for identifying 
a transient internal characteristic value in a transient state of the functional model by using the steady internal char- 
acteristic value (claim 1). 

[0032] Thus, in the characteristic value identification method of the present invention, the functional model of the 
part is prepared by the energy rules called the potential quantity and the flow quantity, and the transient identification 
is performed after the steady identification of its internal characteristic value, thereby enabling accurate identification 
of the characteristic values for all the parts and modeling of an integrated product. 

[0033] The above-mentioned second process may include a first step for determining an internal characteristic value 
of at least one steady test model from the steady functional model, a second step for collecting steady test data by 
performing a test corresponding to the steady test model, and a third step for identifying a steady internal characteristic 
value of the internal characteristic value based on the steady test data (claim 2). 

[0034] The above-mentioned first step may determine the internal characteristic value from a government equation 
in the steady state of the functional model (claim 3). 

[0035] Also, the above-mentioned third step may convert the government equation into a recurrence equation to 
determine the steady internal characteristic value from a recurrence coefficient of the recurrence equation (claim 4), 
and may divide the steady internal characteristic value into a known factor and an unknown factor to identify the steady 
internal characteristic value of the unknown factor (claim 5). 

[0036] Furthermore, the above-mentioned third process may include a first step for determining an internal charac- 
teristic value of at least one transient test model in a transient state of the functional model, a second step for collecting 
transient test data by performing a test corresponding to the transient test model, a third step for applying the steady 
internal characteristic value to the internal characteristic value of the transient test model to generate transient phe- 
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nomenon reproduction data, and a fourth step for correcting the transient phenomenon reproduction data based on 
an error between the transient phenomenon reproduction data and the transient test data, thereby identifying a transient 
internal characteristic value (claim 6). 

[0037] When the error does not lie within an allowable range, the above-mentioned fourth step may repeatedly correct 
a predetermined transient internal characteristic value within the transient phenomenon reproduction data until the 
error lies within the allowable range, and determine the transient internal characteristic value to be identified when the 
error lies within the allowable range (claim 7). 

[0038] Also, the above-mentioned fourth step may preliminarily calculate a variance deviation, as a time history 
sensitivity, to an initial value at a time when each transient internal characteristic value is increased or decreased at a 
fixed ratio, and select a transient internal characteristic value having a maximum sensitivity within the time history 
sensitivity as the predetermined transient internal characteristic value (claim 8), or may select a transient internal 
characteristic value having the time history sensitivity similar to the error as the predetermined transient internal char- 
acteristic value (claim 9). 

[0039] It is to be noted that the fourth step may simultaneously select a plurality of transient internal characteristic 
values having different maximum sensitivity times as the predetermined transient internal characteristic value (claim 
10). 

[0040] Also, a characteristic value identification apparatus according to the present invention comprises block re- 
placement means for a functional model of a part, test reproduction means for reproducing at least one steady test 
model in a steady state of the functional model and at least one transient test model in a transient state, testing means 
of the part for performing a steady test and a transient test respectively corresponding to the steady test model and 
the transient test model, measurement means for collecting steady test data and transient test data at a time when a 
steady test and a transient test of the part are performed by the testing means, and calculating means for identifying 
a steady internal characteristic value of the steady test model by using the steady test data, for applying the steady 
internal characteristic value to the transient test model to generate transient phenomenon reproduction data, and for 
correcting the transient phenomenon reproduction data based on an error between the transient phenomenon repro- 
duction data and the transient test data, thereby identifying a transient internal characteristic value (claim 11). 
[0041] Thus, it becomes possible to quickly identify the functional model of the same kind. 
[0042] When the error does not lie within an allowable range, the above-mentioned calculating means may repeatedly 
correct a predetermined transient internal characteristic value within the transient phenomenon reproduction data until 
the error lies within the allowable range, and determine the transient internal characteristic value to be identified when 
the error lies within the allowable range (claim 12). 

[0043] Also, the above-mentioned calculating means may preliminarily calculate a variance deviation, as a time 
history sensitivity, to an in itial value at a time when each transient internal characteristic value is increased or decreased 
at a fixed ratio, and select a transient internal characteristic value having a maximum sensitivity within the time history 
sensitivity as the predetermined transient internal characteristic value (claim 13), or may select a transient internal 
characteristic value having the time history sensitivity similar to the error as the predetermined transient internal char- 
acteristic value (claim 14). 

[0044] It is to be noted that the calculating means may simultaneously select a plurality of transient internal charac- 
teristic values having a different maximum sensitivity time as the predetermined transient internal characteristic value 
(claim 15). 

[0045] Furthermore, a virtual testing system according to the present invention incorporates the functional model, 
as a virtual prototype, having an internal characteristic value identified by a characteristic value identification apparatus 
claimed in claim 11 and comprises condition assigning means for assigning a driving operation condition and an en- 
vironment condition to the characteristic value identification apparatus, observation means for observing reproduction 
data obtained by the virtual prototype when the driving operation condition and the environment condition are assigned, 
and evaluation means for evaluating an observation result of the observation means (claim 16). 
[0046] Namely, in the identification of the internal characteristic value by the present invention, an evaluation of a 
function, a performance, and a characteristic of a part as well as a product in an actual machine test can be regarded 
as a virtual test performed on a computer, and the actual machine test performed in a development process of designing, 
prototyping, and testing in the prior art can be omitted, so that it becomes possible to shorten the period and to reduce 
the development cost by the virtual test. 

[0047] The above-mentioned virtual testing system may further comprise another measurement means for measuring 
actual machine test data at a time when the driving operation condition and the environment condition are provided to 
an actual machine which forms a subject of the virtual prototype, and re-identification means of the virtual prototype, 
wherein the evaluation means may compare an output of the measurement unit and the observation result, and make 
the re-identification means re-identify the virtual prototype according to the comparison result (claim 1 7). 
[0048] A fixed virtual prototype may be incorporated into a part of a drive system and a load system connected to 
the part as the above-mentioned virtual prototype, the testing means may perform a test corresponding to the fixed 
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virtual prototype, and the evaluation means at this time may make the re-identification means perform a re-identification 
according to the comparison result (claim 18). 

Brief Description of Drawings 

[0049] 

Fig.1 is a block diagram showing a basic concept of a characteristic value identification method and an apparatus 
therefor according to the present invention; 

Fig.2 is a block diagram showing a basic system of a motor model with a potential quantity and a flow quantity of 
energy being mutually related; 

Fig.3 is a diagram showing a general electric circuit of a motor; 

Fig.4 is a diagram respectively showing a functional model and symbols of a motor used in a characteristic value 
identification method and an apparatus therefor according to the present invention; 
Fig.5 is a diagram showing a general electric circuit of a battery; 
Fig.6 is a diagram showing a functional model of a battery; 

Fig.7 is a block diagram showing a process for preparing a steady test model group from a motor steady functional 
model relating to a characteristic value identification method and an apparatus therefor according to the present 
invention; 

Fig.8 is a diagram showing a steady functional model of a motor; 

Fig.9 is a diagram showing a motor model at the time of a steady load test; 

Fig.1 0 is a diagram showing a motor model at the time of a load test; 

Fig. 11 is a diagram showing a motor model at the time of a lock test; 

Fig. 12 is a diagram showing a motor model at the time of an induced voltage test; 

Fig. 13 is a diagram showing a motor model at the time of a braking test; 

Fig. 1 4 is a block diagram showing an identification process of a steady functional model from a steady test model 
group in a characteristic value identification method and an apparatus therefor according to the present invention; 
Fig.1 5 is a graph showing a motor characteristic of a steady state; 
Fig. 16 is a diagram showing an external shape of a motor; 

Fig. 1 7 is a graph showing a steady internal characteristic value of a motor with an ordinate and an abscissa in Fig. 
15 being exchanged; 

Fig.1 8 is a graph showing a steady internal characteristic value of a motor; 

Fig.1 9 is a block diagram showing a transient identification process of a functional model in a characteristic value 

identification method and an apparatus therefor according to the present invention; 

Fig.20 is a diagram showing a structure of an electrically driven arm mechanism; 

Fig.21 is a diagram showing a model at the time of an actual load transient test of a motor; 

Fig.22 is a diagram showing a transient test model at the time when a motor axis is fixed; 

Fig.23 is a diagram showing a fixed torque transient test model of a motor; 

Fig.24 is a diagram showing a model at the time of an inertial rotation test of a motor; 

Fig.25 is a diagram showing a motor operation switch circuit; 

Fig.26 is a diagram showing a functional model of a motor operation switch; 

Fig.27 is a diagram showing an electric circuit of a switching method; 

Fig.28 is a diagram showing a functional model of an electric circuit of a switching method; 

Fig.29 is a diagram showing an electric circuit of a fixed current control; 

Fig.30 is a diagram showing a functional model of a fixed current control circuit; 

Fig.31 is a graph showing a trial simulation result before performing a transient identification in a characteristic 
value identification method and an apparatus therefor according to the present invention; 
Fig.32 is a flow chart illustrating an operation of a transient identification in a characteristic value identification 
method and an apparatus therefor according to the present invention; 

Fig.33 is a graph respectively showing a sensitivity of a characteristic value in a characteristic value identification 
method and an apparatus therefor according to the present invention; 

Fig.34 is a graph respectively showing a relationship between a current deviation and a characteristic value time 
history sensitivity in a characteristic value identification method and an apparatus therefor according to the present 
invention; 

Fig.35 is a graph respectively showing a simulation result after performing an identification in a characteristic value 
identification method and an apparatus therefor according to the present invention; 

Fig.36 is a graph respectively showing a simulation result of a start and a stop in a characteristic value identification 
method and an apparatus therefor according to the present invention; 
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Fig.37 is a graph respectively showing a simulation result of a regular/reversed rotation and brake in a characteristic 
value identification method and an apparatus therefor according to the present invention; 
Fig.38 is a block diagram illustrating a characteristic value identification apparatus according to the present inven- 
tion by taking a motor as an example; 
5 Fig. 39 is a block diagram showing the relationship between a product development and a model; 

Fig.40 is a block diagram showing a concept of a virtual test according to the present invention; 
Fig.41 is a diagram showing a model at the time when a virtual testing system according to the present invention 
is applied to a car; 

Fig.42 is a block diagram showing an embodiment of a virtual testing system according to the present invention; 
10 Fig.43 is a block diagram showing an example applying a virtual testing system according to the present invention 

to an actual machine test of a transmission; 

Fig.44 is a model diagram showing a virtual testing system of a drive system, a transmission, and a load system; and 
Fig.45 is a block diagram showing an example applying a virtual testing system according to the present invention 
to an actual machine test of a power train. 

15 

Description of Reference Numerals 
[0050] 

1 motor model (block replacement structure) 

2 test reproduction model 

3 tested motor 

4 drive circuit unit 

5 load generator 

6 electric load 
8, 34 measurement unit 

9 calculation unit 

1 0 virtual test model 

11 virtual prototype 

12 test standard model 

13 engine model (motor model) 

14 body model (generator model) 

21 drive model (battery operation model) 

22 load model (load generator model) 
23, 37, 40 observation model, observation system 

24 characteristic value update model 

25 model replacement unit 
30 virtual test portion 

31 , 38, 39 drive operation model, operation system 

32, 41 environment condition model, environment system 
33 actual machine test portion 

35 evaluation model 

36 control system input/output unit 

50 tested transmission 

51 tested power train (P/T) 

71 drive system control unit 

72 load system control unit 

[0051] Throughout the figures, like reference numerals indicate like or corresponding components. 

so 

Best Mode for Carrying Out the Invention 

[0052] Hereinafter, an embodiment of a characteristic value (model) identification method and an apparatus accord- 
ing to the present invention schematically shown in Fig.1 will be described. 

55 

1 . Functional model in steady state 

[0053] In order to perform the steady identification process 2 shown in Fig.1 , the steady functional model conversion 
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is required as shown in the process 21 . In this steady functional model conversion, it is necessary as will be described 
hereafter to convert the functional model into the steady functional model for reproducing the steady state of the part, 
and to further convert the steady functional model into the steady test model according to steady tests of the part 
(process 210). 

s [0054] Since the steady identification utilizes the steady test data, the steady test model can make the steady func- 
tional model of the identification subject accord with the tested object to which the actual machine test is performed, 
by modeling the tested object according to the condition and the contents of the steady test for generating the test 
data. Fig.7 shows a flow of a series of model conversions. 

10 1 .1 Conversion from function model into steady functional model 



[0055] Since there is no transient variation in a steady state of a motor, neither the inductance L M nor the moment 
of inertia J M shown in Fig. 4A influences the steady internal characteristic vaiue. Also, since the insulation resistance 
R c in Fig. 4A has an extremely large value, the influence to the motor current l M can be neglected. From the above- 
15 mentioned points, the motor model in the steady state where L M , J M> R c of Fig.4(1 ) are eliminated turns into a functional 
model shown in Fig.8. 

[0056] The left side of Fig.8 shows a battery model of Fig. 6 in which the operation switch of the electric system is 
omitted, and the right side thereof shows a motor model in a steady state. From Fig.8, the motor model in the steady 
state can be expressed by the following government equation: 
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[0057] For another method of deriving the steady government equation, the government equation in the steady state 
obtained by eliminating L M ,J M , from the government Eq.(3) of the functional model in Fig.4(1) can be obtained by 
the following processes: 

® The equation obtained by eliminating the differential state values of x M and x M from Eq.(3) with L M and J M 
being assumed to be 0 is given as follows: 
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@ The equation obtained by eliminating the internal state values of x M and X|_ from equation (9) is given as follows: 
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"D" in equation (10) is assumed to be the following equation: 

d = d m( r m+ r c)+8 m M m 2 Eq.(11) 
® The equation obtained by eliminating the insulation resistance by making R^oo from Eq.(1 0) is given as follows: 
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[0058] The government equation (1 2) in the steady state derived from the government Eq.(3) of the functional model 
is equal to the government equation (8) derived from the functional model in the steady state in Fig.8. 
[0059] Accordingly, there are two methods for obtaining the steady government equation reproducing the steady 
state of the part; one by converting the functional model reproducing the transient state of the part into the steady 
functional model, and the other by deriving the steady government equation from the government equation of the 
functional model. This indicates that the functional model in the steady state can be reversely depicted from the derived 
steady government equation, and that there is a reversibility between the steady functional model and the steady 
government equation. 



1 . 2 Conversion from steady functional model into test model 



[0060] A steady state of an actuality by the steady functional model will be considered, and the test model will be 
examined. The test model in the steady state models the test of the actuality which provides a specified condition to 
an input/output system of a tested object, from the steady functional model and its government equation. The steady 
internal characteristic value is identified by relating the test data obtained when the specified condition is assigned to 
the tested object, with the government equation. 

[0061] The relationship between the specified condition in which the steady internal characteristic value is assigned 
to the motor model and the test is as follows: 

© In case a fixed torque is provided to the motor with the power supply being connected, the model reproduces 
the relationship between the load applied to the motor and the steady internal characteristic value, so that the 
relationship between the motor current, the angular acceleration, and the load torque can be observed (steady 
load test: Fig.9). 

© The case where the input torque T M is assumed to be 0 corresponds to a motor rotation in an idle state, so 
that an idle angular velocity and an idle current are observed (idle test: Fig.tO). 

® The case where an output angular velocity is assumed to be 0 corresponds to a lock state or a state at a 
starting moment when the motor rotation is compulsorily stopped, so that a lock current and a lock torque can be 
observed (motor lock test: Fig.11). 

@ In case the input current l M is assumed to be 0 and the input torque T M is provided, the model which reproduces 
a generator function of the motor is obtained, so that the induced voltage and the angular acceleration can be 
observed (induced voltage test: Fig.12). 
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(D In case a resistor is connected to the input/output state value of the electric system to provide a fixed torque, 
the model which reproduces a regenerative braking force of the motor can be obtained, so that the relationship 
between the motor current, the angular acceleration, and a braking torque can be observed (steady braking test: 
Fig. 13). 

[0062] These test models are different from each other depending on how to provide the specified conditions. How- 
ever, in an actual modeling, the values of the steady characteristic such as a steady functional model or a steady load 
test model and the input state are made 0 or an infinite value (large value within a range where the value does not 
influence a reproduction result), so that each test model can be equivalent^ realized. It is necessary to perform at least 
one of such steady tests. 

[0063] In view of the above-mentioned description, the motor test model performed from the steady functional model 
in Fig.8 under the specified conditions, is given as follows: 

1 . 2. 1 Motor steady load test 

[0064] Fig.9 shows a steady load test which observes the performance and the characteristic by providing the load 
torque to the motor. 

[0065] In Fig.9, the input torque T M is provided as a drive torque by the additional load from the steady functional 
model in Fig.8, and the output angular velocity is represented as an observation quantity. Accordingly, the govern- 
ment Eq.(8) of the motor is transformed into the following equation: 



M 



M 



M 



•LA 



K M + — E MB' 



M 

M M 5 M 



(T M + T MF ) 



...Eq. (13) 



[0066] The government equation of the steady load test obtained by mutually substituting the government Eq.(4) of 
the battery for equation (13) is given as follows: 



OX 



M 



'M 



'M 



1 + 

M.,S, 



'M 



R B +R M 



M M 2 V 



.Eq.(14) 



[0067] Equation (1 4) is a government equation reproducing a rated load test of the motor, where all of the equations 
assume the output state value reproducing the state value of the steady state. 
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1 . 2. 2 Motor idle test 

[0068] In the motor idle state, the input torque T M of the steady load test model in Fig. 9 becomes 0. The idle test 
model where the input torque T M is eliminated from Fig.9 can be shown by Fig.1 0. 

[0069J The government equation of the motor in Fig. 10 turns into the following equation obtained by eliminating the 
load torque T M from the government Eq.(13) of the steady load test: 
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...Eq.(15) 



[0070] The government equation of the idle test obtained by mutually substituting the government Eq.(4) of the battery 
for equation (15) is given as follows: 



1 m mCVE mb )-(R b +R m )T mp 
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M 
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1 + 

, M M 5 M T 



MF 



'M 



...Eq.(16) 



1 .2.3 Motor lock test 



[0071] At the starting moment of the motor, the angular velocity becomes 0, where it can be regarded as a motor 
lock state in which the motor rotation is compulsorily stopped. Accordingly, the functional model of the lock test repro- 
ducing the motor lock can be shown by Fig. 11 where the output angular velocity of the steady functional model in 
Fig. 8 is supposed to be 0. It is to be noted that since a brush voltage drop E MB generated in the stop state of the motor 
rotation does not generate a fixed voltage drop but only assumes a minute voltage drop by a brush resistance, it is 
neglected and omitted. 

[0072] In Fig. 11 , each steady internal characteristic value relating to the output angular velocity is eliminated 
from the steady functional model in Fig.8, so that the input torque T M is shown by the observation quantity of the lock 
torque. The government equation of Fig. 11 is given as follows: 



55 



[I M ] ...Eq.(17) 
[0073] The government equation of the lock test obtained by substituting the government Eq.(4) of the battery for 
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equation (1 7) to be integrated is given as follows: 




M 




[1] 



Eq.(18) 



M 



R B +R M 

E 0 



R B +R 1 



1 .2.4 Motor Induced voltage test 

[0074] Fig. 12 shows a model of an Induced voltage test which provides the torque to the motor to observe a counter 
electromotive force V M induced by the armature and the angular velocity c^. 

[0075] In Fig.12, the steady internal characteristic value relating to the input current l M is eliminated from the steady 
functional model in Fig.8, so that the output angular velocity is shown by the observation quantity. The government 
equation of Fig.12 is given as follows: 



1 .2.5 Motor steady braking test 

[0076] Fig.1 3 shows a steady model of a regenerative braking which provides the torque to the motor and reversely 
returns the counter electromotive force V M induced by the armature to the motor through a braking resistance R L to 
provide a braking force. 

[0077] In Fig.13, the battery is eliminated from the steady functional model in Fig.8, the braking resistance R L is 
connected, and the input torque T M of the motor is provided as a drive torque. In addition, the output angular velocity 
com is shown as the observation quantity. The government equation of the motor in Fig.13 is given as follows: 




[i] 



...Eq.(19) 
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MmEmb + CRl+RmH'VT'mf) 1 



0), 



'M 



M 



'M 




[1] 



...Eq.(20) 



M 



M M 26 M 



[0078] The above-mentioned test models can be equivalently modeled by setting the following values to the char- 
acteristic or the input state value of the steady load test model shown in Fig.9: 

© The idle test model becomes equivalent by rendering the input torque T M = 0. 

@ The lock test model becomes equivalent by rendering a viscous damping coefficient and the brush voltage 
drop respectively D M = E MB = 0. 

® The induced voltage test model becomes equivalent by rendering the motor current and the brush voltage drop 
respectively l M = 0,E MB = 0. 

@ The braking test model becomes equivalent by rendering the internal electromotive force of the battery E 0 =0, 
and by rendering the internal resistance Ffe the braking resistance R L . 

2. Identification of steady state 

[0079] Fig. 14 shows an identification process 2 of the steady internal characteristic value included by the steady 
functional model having been described referring to Fig.1 . Since the steady internal characteristic value known before- 
hand from a design value, the characteristic value of a similar article, the measured value of the characteristic, and 
the like can be determined as a fixed value in the steady identification process (process 26), remaining unfixed values 
can be identified (process 27). 

[0080] As for an analysis technique applied to this, the government equation of the steady test model is converted 
into a statistical model such as a canonical correlation analysis or a multiple fold recurrence analysis (simple recurrence 
analysis), determines a deviated recurrence coefficient from the steady test data, and derives the steady internal char- 
acteristic value from the relationsh ip between the deviated recurrence coefficient and the government equation (process 
28). Also, as for the steady characteristic formed by many factors, the factors are analyzed, and divided into factors 
of which quantification is possible and factors of which quantification is difficult by modeling the forming process of the 
steady internal characteristic value with a theoretical equation, an experimental equation, and the like, thereby identi- 
fying the factors of which quantification is difficult (process 29). 

2. 1 Test plan of steady identification 

[0081] In the steady functional model shown in Fig.8, there are six kinds of the steady internal characteristic values 
m m> r m- d m» &m> t mf» and e mb wn ' cn for7TI tne subject of the steady identification, in order to determine these steady 
internal characteristic values, it is necessary to firstly clarify the function of the steady model in Fig.8 and the charac- 
teristic of the steady internal characteristic values, and to clarify the relationship with the test method. 

2.1.1 Verification of known/measurable steady internal characteristic value (process 26) 

[0082] Firstly, known characteristic values and the characteristic values which can be solely measured are clarified 
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in the steady internal characteristic values included In the steady functional model. Among the characteristic values 
which govern the steady state of a general product or part, there are a large number of characteristic values which can 
be clarified without performing the steady test such as a characteristic value clarified in the literature or the like, a 
characteristic value which can directly utilize the design value, and a characteristic value which can be directly meas- 
ured. Accordingly, it is important to determine the characteristic values by a preliminary investigation or measurement 
in case of the steady identification. 

[0083] Reviewing the motor used here, such characteristic values to be determined are a brush voltage drop E MB , 
an insulation resistance Rq omitted at the time of modeling the steady functional model in Fig.8, a motor winding 
resistance R M , and the like. These steady characteristics are given as follows: 

® Since being a voltage drop due to a graphite brush, the brush voltage drop E MB assumes constant regardless 
of the current l M during the motor rotation, and assumes a minute DC resistance in the stop state. Accordingly, 
0.7 [V] shown in the literature is applied during the motor rotation, and the direct current in the stop state can be 
neglected. 

(D The insulation resistance R c assumes the value [MQ] which does not influence the result. 

® The winding resistance R M of the motor can be determined by the ohmmeter or by calculating the motor voltage 

and current in the motor lock state shown in Fig. 11 . 

2.1.2 Steady functional model and statistical model (process 28) 

[0084] There is a multivariate analysis such as a multiple fold recurrence analysis and a simple recurrence analysis 
in the analysis technique frequently used for the product development. The recurrence equation of the multiple fold 
recurrence model applied thereto is given as follows: 

y = b 0+ b 1 x 1 +b 2 x 2 ....b n x+e rr Eq.(21) 

[0085] Applying the multiple fold recurrence equation shown in equation (21 ) to the steady functional model including 
a plurality of unknown steady internal characteristic values will be examined. As for the method of application, depend- 
ent variables and independent variables in the steady government Eq.(8) are made correspond to an objective variable 
y and predictor variables x r x n , to be determined from the steady test data. 

[0086] Also, deviated recurrence coefficients b 0 -b n are values combining known and unknown steady internal char- 
acteristic values. Accordingly, unknown steady internal characteristic values can be derived by relating (associating) 
the deviated recurrence coefficients to (with) the steady internal characteristic values. In Eq.(21 ), e^ is a residual. Also, 
when the multiple fold recurrence equation has a plurality of different objective variables, these are expressed by 
simultaneous equations so that the canonical correlation analysis is applied. As a matter of course, the test data ob- 
tained by a plurality of different test methods, the deviated recurrence coefficient, and the steady internal characteristic 
value of the recurrence model can be related to each other. 

[0087] In the multiple fold recurrence analysis of the motor, the steady internal characteristic value can be clarified 
by the following process from the data of the steady load test shown in Fig.9: 

[0088] Firstly, as for the multiple fold recurrence model, the following multiple fold recurrence equations can be de- 
rived from the government Eq.(13) of the steady load test model in Fig.9: 

T M = t mf- d m w m +m m'm Eq.(22) 

T m = -(T Mp + — 1 +I 7T-r- + M M ll M ...Eq.(23) 

v M M 6 M > mm ' 

[0089] It is to be noted that equation (22) corresponds to the 1st row of Eq.(13), and equation (23) corresponds to 
the 2nd row of Eq.(1 3). Since the internal resistance R B of the battery is very small, and the internal electromotive force 
E 0 of the battery and a motor output voltage V M are considered to be approximately equal, E 0 in Eq.(23) is substituted 
for the dependent variable V M (motor output voltage) of Eq.(1 3). 

[0090] The recurrence equation with the torque T M of these equations being made an objective variable and the 
motor current l M and the angular velocity being made predictor variables is given as follows (process 281 ), where 



15 



EP 1 288 759 A1 



the upper equation of the following equations (24) forms a multiple fold recurrence equation with two predictor variables, 
and the lower equation forms a simple recurrence equation with a single predictor variable: 



m a B 10 +B n I M + B 12 o> M ] _ 

m - b 20 *b 21 i m | ~*«W> 



T = 
T, 



[0091 J Also, constants (Intercepts) B 10 , B^, and deviated recurrence coefficients B 11f B 12 ,B 21 can be expressed by 
the following equations: 



B,o 


= -T 

MF 


B„ 


= -Dm 


B 12 


= M M 


B 20 


D E 0 _E MB 
M M M 6 M 


B « 


= M M *-^-^ 

M M 6 M 



-T, 



MF 



...Eq.(25) 



[0092] The steady internal characteristic values derived from the deviated recurrence coefficients shown in equations 
(25) are given as follows (process 282): 

[0093] Firstly, T MF , D M , M M correspond to the respective deviated recurrence coefficients from the upper four rows 
of Eqs. (25). 



MF 



-B 



10 



D M = - B ll 



B, 



..Eq.(26) 



*M "12 

[0094] Then, ^ can be derived by substituting equation (26) for four rows of Eqs. (25) as follows: 

^'(Baa-BtoJBia <E ° ' Emb) 



Eq.(27) 



[0095] It is to be noted that the induced voltage E 0 of the battery and the brush voltage drop E MB in equation (27) 
are assumed to be known steady internal characteristic values. 

[0096] Moreover, R M can be derived from the following equation by substituting Eqs. (26) and (27) for five rows of 
Eqs. (25): 



R " = B^ (E o- E mb) 



Eq.(28) 



[0097] While in the above-mentioned example, ail of the steady internal characteristic values have been determined 
from the government Eq.(8) by the recurrence analysis, in many cases of the actual steady identification, the steady 
internal characteristic values are derived by combin ing a plurality of test methods. For example, in the above-mentioned 
example, as for the 3rd row of Eqs.(26) deriving the motor constant M M and equation (28) determining the winding 
resistance R M , both M M and R M can be determined from Eq.(1 7) of the lock test model and the test data as well. 
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2.1 ,3 Factor analysis of steady characteristic (process 29) 

[0098] The steady internal characteristic value of the steady functional model includes many internal factors forming 
this characteristic. Accordingly, it is necessary to clarify the factors composing each characteristic value by performing 
a factor analysis, and to formulate the result. Also, it is necessary to divide the factors into known factors and unknown 
factors based on the result, and to perform a steady identification for clarifying an influence degree of the unknown or 
uncertain factors by the model or the test data. 

[0099J To the steady identification method, a sensitivity analysis or the like is effective by which the relationship 
between the variance of the government factor and the steady internal characteristic value is taken as a sensitivity. 
Particularly, in case the steady internal characteristic value has a non-linear characteristic influenced by another state 
value, the sensitivity analysis is important. 

[0100] If the motor constant is considered, for example, P, "a", and Z are known values in Eq.(1 ) since they are design 
constants determined by the motor structure. However, it is difficult to determine all of the magnetic fluxes $ per pole 
since it is distributed in the gap between the field and the armature. Accordingly, to identify the motor constant means 
to identify the magnetic flux $ by the test data. All of the magnetic fluxes <j> per pole determined from the identified motor 
constant M M are given as follows: 

O = M M ?gx10 8 Eq.(29) 

[0101] Moreover, as a brief example, the winding resistance R M of the motor is given by the following equation 
determined by a line length "I", a sectional area "S", a specific resistance "p", and a temperature coefficient V, being 
governed by a temperature T: 

r m = P^0+<*t) Eq.(30) 

[0102] From equation (30), I, S, p, and a assume the design values determined by the structure and the stuff, while 
the temperature V is governed by the temperature of the test environment and the temperature rise by self-heating 
of the motor. Accordingly in the test for performing the identification, the temperature of the motor forms an important 
factor which governs the characteristic value. Also, upon the examination of the high temperature environment and 
low temperature environment, it is required to perform a temperature correction, by Eq.(30), to R M in the government 
Eq.(3) and the steady government equations of equation (8) and subsequent equations. 

2.1 .4 Motor steady characteristic 

[0103] The steady internal characteristic value of the motor is represented by a torque-current characteristic repre- 
senting the relationship between the output torque T M [Nm] and the motor current l M [A], and by a torque-velocity 
characteristic representing the relationship with the angular velocity co M [rad/sec], based on the result measured in the 
above-mentioned steady load test (Fig.9). Firstly, the torque-velocity characteristic can be expressed by the following 
equation obtained by substituting the motor current l M of Eqs.(22) and (23), the latter itself representing the torque- 
current characteristic: 



M M ' V R M 7 

[0104] If the torque-current characteristic and the torque- velocity characteristic of the motor are available from a 
catalog, technical materials, or the like, the steady internal characteristic value may be determined from the charac- 
teristic diagram. This motor characteristic is shown in Fig. 15. 

[0105] The relationship between the two recurrence lines in Fig. 1 5 is expressed by Eq. (22). In Fig. 15, the intersection 
of X axis and the torque-angular velocity characteristic shows the motor lock state, and the intersection of Y axis and 
the torque-angular velocity characteristic shows the idle state. Accordingly, as for the motor constant M M and the 
winding resistance R M , their approximate values can be determined by substituting the motor current l M , the torque 
T M , and the voltage V M in the lock state for Eq.(1 7). Also, by substituting the result for Eq.(22), an approximate value 
of the viscous resistance coefficient 6^, at the time of idling can be determined, whereas the friction torque T MF is 
neglected. 
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2.2 Motor steady test and steady identification result 

[0106] The characteristic value in the steady state is identified from the motor functional model shown in Fig.4(1 ). It 
is to be noted that the motor used in this identification is a small DC motor of a separately excited type widely used for 
an electric actuator of a car, and a permanent magnet is used for the field. 

[0107J The steady test for the steady identification has been performed to the idle test, the lock test, and the steady 
load test. The idle test and the lock test have been performed as a test under the specified conditions of the motor. 
The former has been performed with the output axis of the motor being released in the state of the input torque T M = 

0, and the latter has been performed with the output axis being fixed in the state of the output angular velocity w^, = 0. 
[0108] Moreover, in the steady load test, a string hung with a weight has been winded onto a drum attached to the 
motor output axis, whereby the load torque has been calculated from the mass of the weight and the radius of the 
drum. Also, the test conditions are shown in the following Table 3 together with the result of the steady identification. 
[0109] Fig.1 7 shows a result of the steady test performed to the motor of Fig.1 6 and the recurrence equations, where 
X and Y axes of Fig. 15 are exchanged. 

[0110] Fig.17 is a diagram showing the measured values of the torque-current characteristic representing the rela- 
tionship between the output torque T M [Nm] in the steady state and the motor current l M [A], and representing the 
relationship with the angular velocity co M [rad/sec], and its recurrence equation. The characteristic value on the X axis 
represents the motor idle state (T M = 0) as mentioned above, and the intersection between the torque-velocity char- 
acteristic and the Y axis shows the lock state or the starting moment of the motor with the angular velocity (0^ = 0. The 
points shown by variable names in Fig.1 7 show a basic characteristic and a rating characteristic (dotted line) in Table 

1 . Also, the black circle in Fig.1 7 shows the measured values. 

[0111] Moreover, the basic performance at the time of starting and idling obtained in Fig.17 is shown in Table 1, and 
a catalog specification published by a motor maker is shown in Table 2 as reference data. 



TABLE 1 





CHARACTERISTIC NAME 


SYMBOL 


UNIT 


CHARACTERISTIC VALUE 


MOTOR 


IDLE CURRENT 


'm_o 


[A] 


0.1 




IDLE ANGULAR VELOCITY 




[rad/sec] 


1410.0 




STARTING TORQUE 


T M_L 


[Nm] 


28.0X10- 3 




STARTING CURRENT 


! M_L 


[A] 


4.0 


TEST POWER SUPPLY 


INTERNAL 

ELECTROMOTIVE FORCE 




[V] 


12.0 




POWER SUPPLY INTERNAL 
RESISTANCE 


R B 


[A] 


1.0X10- 3 



TABLE 2 



CHARACTERISTIC NAME 


SYMBOL 


UNIT 


CHARACTERISTIC VALUE 


CONDITION 


RATED CURRENT 


'm_s 


[A] 


0.57 


MAX 


0.9 


RATED ANGULAR 


<*>M_S 


[rad/sec] 


1210.0 


TOLERANCE 


±160.0 


VELOCITY 












RATED LOAD 


T M S 


[Nm] 


4.0 x 10-3 






RATED VOLTAGE 




[V] 


12.0 







[01 12] Table 3 shows a result of the steady identification by the test data of Fig.1 7. 
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TABLE 3 





CHARACTERISTIC 
NAME 


SYMBOL 


UNIT 


CHARACTERISTIC 
VALUE 


MOTOR STATE 


VISCOUS RESISTANCE 




[Nm/(rad/sec)] 


0.65 x10" 6 




KINETIC FRICTION 
TORQUE 


! T MF 


[Nm] 


0.1 x10* 6 




WINDING RESISTANCE 
BRUSH VOLTAGE 


Rm 


[H] 
[V] 


3.0 
0.7 




INSULATION 
RESISTANCE 




[Mfi] 


10.0 




MOTOR CONSTANT 


M M 


[Nm/A] 


7.25x10-3 




VELOCITY 

COEFFICIENT 

CORRECTION 


5m 


M 


1.0 


TEST CONDITION 
(KNOWN) 


TEST VOLTAGE 

INTERNAL 
RESISTANCE 


E 0 
Rb 


V 


12.0 

50.0X10*3 




AMBIENT 
TEMPERATURE 


t F 


tdeg] 


25.0 



[0113J Table 3 is a result assembled according to the process of the steady identification examined so far. Among 
the identified steady internal characteristic values, E BR has been assumed to be a prescribed value determined from 
the literature. Also, Rc is a measurement result of the insulation testing set. The other steady internal characteristic 
values M M , R M , D M 5 M , and T MF are determined from the multiple fold recurrence analysis and the simple recurrence 
analysis by the above-mentioned Eqs.(24)-(27). Also, since the internal resistance R B of the battery is so small that 
the resistance value of the current measurement distributary unit has been applied to the internal resistance R B of the 
battery. 

[0114] It is to be noted that in order to determine five coefficients B 10 , B 20 , B l1t B 12 , and B 21 , five simultaneous 
equations are required. T M , l M> and are vectors where equal to or more than five values determined by the char- 
acteristic of Fig. 1 7 exist. 

[0115] Fig.1 8 shows a simulation result of the thus identified steady internal characteristic values by using the steady 
government equation obtained by substituting Eq.(4) of the battery model for Eq.(8) of the motor steady state. 
[0116] Thick solid lines in Fig. 18 show results of the steady state of the angular velocity and the motor current 
l M , while thin dotted lines show simulation results performed in Eq.(5) by using the transient internal characteristic 
value before identification described later. It is seen from Fig.1 8 that the results of the steady state coincides with each 
other in the steady region where the transient state becomes stable. 

3. Identification of transient state (process 32) 

[0117] Since the steady internal characteristic values of the motor shown in the above-mentioned Table 3 are clarified, 
the moment of inertia J M and the inductance L M as transient internal characteristic values governing the transient 
function of the motor are identified next with the measured value of the motor current l M , by the transient identification 
process 32 shown in Fig.1 9. 

[0118] As an identification method, approximate values are firstly provided to the transient internal characteristic 
values to perform a simulation so that initial values are determined (process 36), then the time history sensitivities of 
the transient internal characteristic values for the motor current l M are determined (process 38), and the deviation 
(process 37) between the measured value (process 35) of the motor current l M and the simulation value assumes the 
least for the sensitivity indexed. 

[01 19] Moreover, the viscous resistance coefficient 6 M and the kinetic friction torque F MF already identified are again 
identified in the transient identification for the confirmation. As for the identification, the assembly/disassembly and the 
replacement by a block replacement method can be performed to the model of the operation system and the load 
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connected to the motor source supply system and the rotation system. 
3.1 Conversion from functional model into transient test model (process 33) 

[0120] The test model in the transient state reproduces an actuality test providing a specified condition to the input/ 
output system of the tested object in the same way as the steady identification by the functional model and its govern- 
ment equation. In the transient identification, the test data obtained by providing the specified condition to the tested 
object, and the simulation result by the test model are related, so that the characteristic value group included in the 
government equation is identified. Also, the test model of the transient identification models the drive or the load to be 
incorporated by the following method: 

© Input/output model for driving 

[01 21 ] The input/output model for driving is a model for supplying driving energy such as a power supply and a torq ue 
in order to drive the tested object. This model occasionally has an operation system or a control system for operating 
the supplying method of the driving energy. For example, the system connected to the electric system of the motor 
corresponds to such a system including a battery of a power supply, an operation switch of a motor, a control unit, and 
the like. As a specific example, there is a braking resistance of the regenerative braking mentioned above, or the like. 

(D Input/output model for load 

[0122] The input/output model for the load is a model for accumulating or consuming the energy such as a charac- 
teristic which forms the load of the tested object and an actual load. This model occasionally has the operation system 
or the control system for operating how to provide the load energy. In the rotation system or the mechanical system of 
the motor, there is a method for connecting a stiffness (spring stiffness), the moment of inertia, the viscous resistance 
coefficient, and the like. 

[0123] It is to be noted that in case the characteristic whose value is large is connected, it can be considered that 
the output axis is fixed. Also, as an expediency for providing the fixed load, there is a method by which the current, the 
torque, and the like are provided as an additional load, and the voltage and the angular velocity which become a pair 
are made an observation quantity. For operating a machine load, the load is switched over by simultaneously using 
models of a clutch, a brake, and the like. 

3.1.1 Actual load transient test model of motor (process 34) 

[0124] In the transient test model which reproduces the state of the practical use, the battery model in Fig.6 is united 
with the electric system of the motor model in Fig.4(1 ), and an actual load model or an equivalent load model is further 
united with the rotation system to be converted into a model for an execution closed from the outside. 

(1) Actual load model of electrically driven arm mechanism 

[01 25] Fig.20 shows a load model of an electrically driven arm mechanism having an operation arm whose operation 

range is restricted by 9 max at a rack side of a worm gear, as an example of a practical load. 

[0126] Fig.21 shows a functional model for an execution in which the electrically driven arm mechanism in Fig.20 is 

incorporated in the battery and the motor model. The electrically driven arm mechanism in Fig.21 has a stiffness Kc 

of a stopper, a reduction ratio N G of a reduction unit, and a contact radius l_c of the stopper. Inside the electrically 

driven arm mechanism, there are a presumed rotation angle 9 R of the arm and an internal state 

[0127] Also, the condition judgment shown by a hexagon frame in Fig.21 is a model which judges an arm operation 

range, and operates a switch variable S^. This functional model has an additional load of torque T L as an external load. 

[0128] The government equations of the electrically driven arm mechanism are given as follows: 
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0), 



0 

0 0 

6t ■ N G //X M dt 



0 N G 0 



N G 0 -N 0 T L 



<*>, 



M 



...Eq.(32) 



[0129J In this government equations (32), the upper equation indicates the government equation of the electrically 
driven arm mechanism shown on the right side of Fig.21 , and the lower equation indicates the presumed rotation angle 



6 L of the arm. Also, the model determining the rotation range of the arm is given as follows: 

if («R^m«) then ( s wc = 1 ) else s wc = 0) 



Eq.(33) 



[0130] Equation (33) is for a condition judgment S ws in which ^=0 when the presumed rotation angle 8 R of the 
arm lies within the rotation range G max , and otherwise S wc = 1 . 

[0131] If the above-mentioned Eq.(5) where the battery and motor are integrated, and the government Eq.(32) of 
the electrically driven arm mechanism are integrated, the following equation can be obtained: 
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...Eq.(34) 

[0132] In equation (34), the 1st-3rd rows form the state equations, and the 4th and the following rows form output 
equations for observing the internal transient state value. It is to be noted that Y 0 and l 0 in the equations are an equivalent 
admittance and a current source respectively shown in equations (6) and (7). 



(2) Axis fixed model of motor 



[0133] The mode! in which the output axis of the motor is fixed can be shown by next Fig.22. This model can be 
replaced with the arm mechanism on the right side of the motor actual load transient test model in Fig.21 . 
[01 34] Fig.22 is a model in which the model of the electrically driven arm mechanism in the actual load transient test 
model of Fig.21 is replaced with stiffness Kq which fixes the motor output axis. The motor torque T MF , in Fig.22 repro- 
duces the lock torque in the state where the output axis of the motor is fixed. It is to be noted that the government 
equation is omitted in this case. 
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3.1.1.1 Fixed torque transient test model of motor 



[0135] There is a fixed torque test for providing a fixed torque as a load in the load test performed in the actual 
machine test. This test model unites the battery model in Fig.6 with the motor model in Fig.4(1) to be converted into 
the model for execution closed from the outside. The converted functional model for the execution is shown in Fig.23. 
[0136] In Fig.23, by uniting the battery and the motor, the motor voltage V M and the current l M form the observation 
quantities. Moreover, since the input torque T M of the motor becomes the load torque generated inside the model, the 
input torque T M becomes the additional load T M , and the output angular velocity becomes the observation quantity. 
It is to be noted that in the idle test of the motor, the additional load T M can be made 0. 

[0137] The government equation for execution of the function model therefor shown in Fig.23 is given as follows: 
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...Eq.(35) 



[0138] In equation (35), the 1 st-2nd rows form the state equations, and the 3rd-5th rows form the output equations. 
It is to be noted that the equivalent admittance Y 0 and the current source l 0 are given by the above-mentioned Eqs. 
(6) and (7). 

[0139] Moreover, the process by which Eq.(35) is derived from the above-mentioned Eq.(5) is as follows: 
[01 40] Firstly, the 1 st and 2nd columns including the moment of inertia J M and the inductance L M are eliminated, and 
the dependent variables at the 1st and the 2nd rows are made xm and xj_ Then, the input torque T M of the independent 
variable in Eq.(5) can be substituted for the additional load at the 1st row the 3rd column in Eq.(5). 



3.1.1.2 Inertia rotation test model of motor 



[0141] With the motor, the motor voltage V M , the current l M , the angular velocity co M of the rotation system, and the 
torque T M can be measured. Among them, the electric system particularly has a characteristic that data with high 
accuracy are easily measured. However in the functional model shown in Fig.23, the motor voltage V M nearly equals 
the electromotive force E Q of the battery, so that it can not be used as the data for identification. 
[0142] However, as for the induced voltage V M generated when the motor power supply is turned off, the test data 
for identification can be measured by the inertia! rotation test. The model of this inertia! rotation test is shown in Fig.24. 
[01 43] Fig.24 shows a model of the inertia) rotation test from the moment when the power supply of the motor in the 
idle state is turned off until the motor is stopped. It is to be noted that the insulation resistance R c is omitted because 
the current is minute. The motor model of the inertial rotation test shown in Fig.24 is given as follows: 
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...Eq.(36) 



[0144] In equation (36), the 1 st-2nd rows form the state equations, and the 3rd-5th rows form the output equations. 
It is to be noted that the transient test model of Fig.23 is used for the example of the transient internal characteristic 
value identification performed here, instead of the inertial rotation test model. 

3.1 .1 .3 Actual load transient test model of motor having operation system 

(1) Operation switch 

[0145] The start, the stop, the regular rotation, the reversed rotation of the motor, and the basic operation of the 
electric brake can be performed by a switch circuit shown in Fig.25. The test model in which the modeled switch circuit 
is incorporated into the actual load transient test mode! of Fig.21 will now be examined. 

[0146] The motor operation model can be incorporated between the voltage V M and the current l M which connect 
the battery and the motor in the motor actual load transient test model shown in Fig.21 . By incorporating the operation 
switch, the actual load transient test model in Fig.21 can reproduce the start, the stop, the regular rotation, the reversed 
rotation of the motor, and the drive state of the electric brake by the switch operation. It is to be noted that the simulation 
results of the model will be described after the identification of the transient state. 

[0147] Fig.25 is a circuit diagram of the operation switch incorporated between the battery and the motor. The func- 
tional model of the switch circuit is shown in the next Fig.26. 

[0148] The resistance R L in Fig.26 is the braking resistance for the electric brake. Moreover, neither the contact 
resistance nor the insulation resistance of the switch are considered. 

[0149] The operation switch in Fig.26 is composed of S ws for starting and stopping the motor, for operating the 
regular/reversed rotation, and S WB for selecting necessity/unnecessity of the electric brake. As for an operation signal 
for operating the switches, ON/OFF of the switches are respectively represented by "I" and "0". The functions of the 
switches are as follows: 



© The switch S ws supplies the power to the motor when S ws = 1 , and stops the power supply when S^ s = 0. 
© As for which controls the regular/reversed rotation, a switchover switch variable S wu _ 0 assumes "1 ■ when 
S WL = 0, so that the current l M of the electric circuit in Fig.23 flows in the direction of the arrow to rotate the motor 
regularly. An S WL1 assumes "1 0 when S WL = 1 , so that the current flows reversely to rotate the motor reversely. 
Also, S WL 0 and Sy^j are the switchover switch variables prohibiting that both of 0 and S WL1 simultaneously 
assume the state of "1" or "2". On the functional model, this operation makes Sy^ onlhe current side switch over 
the polarity of the output current l s , so that on the voltage side restores the input voltage V s which has come 
back to the original polarity. 

(D As for the braking by the electric brake, the electric brake is made ON when S WB = 1 and OFF when S WB = 0 
under the condition of the power supply being made OFF, i.e. S ws = 0. The electric brake consumes the induced 
current caused by the motor idling when the power supply is made OFF by the braking resistance R L to generate 
the braking torque. This braking torque is determined by the braking resistance R L , so that the maximum braking 
torque is generated when the resistance value is 0 [CI]. 

[0150] The government equation determined from Fig.26 is given as follows: 
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wS' WE 



...Eq.(37) 



WE 



[0151) S WE of equation (37) is a power supply switch variable for switching over the motor start/stop, and regular/ 
reverse direction of the current flow, which is expressed by the following equation: 



S WE - S WS (S WL _ 0 -S WIO ) 



Eq.(38) 



[0152] The government equation of the power supply system obtained by mutually substituting the voltage V B and 
the current l B of the government equations of the battery and the operation switch in Eqs.(4) and (37) to be integrated 
is given as follows: 



J_ s 



WA 



'WA 



...Eq.(39) 



[0153] S WA in equation (39) is assumed to be a braking switch variable by which ON/OFF of the power supply and 
a necessity/unnecessity of the electric brake are selected, and is expressed by the following equation: 



S WA " S WE 2 + D BS WB^" S Ws) 
OR 

R B 

S WA = S WS + ~ S WB^ "" S Ws) 



.Eq.(40) 



[0154] In equation (40), as for S WE 2 in the upper equation, the switchover switch variable member of the regular/ 
reversed rotation assumes (S^ L _ 0 - S^)* = 1, and S ws 2 assumes either "0" or M 1" from equation (38). Therefore, 
the upper equation can be expressed by "the lower equation arranged by S ws . 

[01 55] Finally, the government equation obtained by mutually substitutin g the motor model of the actual load transient 
test model in Eq.(32) and the voltage V M and the current l M of and the power supply system in Eq.(39) to be integrated 
is expressed by the following equation: 
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[0156] In equation (41), the 1st-3rd rows form the state equations, and the 4th and after rows form the output equa- 
tions. The abstracted equivalent admittance Y 0 and the current source l 0 in Eq.(41) are expressed by the following 
equations: 



0 Rb +r c s wa 



Eq.(42) 



30 



\ 0 - Y 0 E 0 S WE 



Eq.(43) 
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[0157] As an applied example of operating the motor by the electronic circuit incorporated into a control unit or the 
like, it is possible to incorporate the electronic circuit into the actual load transient test model in Fig.21 in the same way 
as the operation switch circuit of the motor in Fig.25. 

[0158] Also, it is possible to replace the operation system of the motor including the above-mentioned switch mech- 
anism by the block replacement method. It is to be noted that while there is another method by which the electronic 
circuit is incorporated by using high circuit technology, the minimum function necessary for driving the motor is modeled 
here. Accordingly, the non-linear characteristic of a semiconductor, the transient response of the electronic circuit and 
the iike are excluded. 



(2) Switching method 



[0159] By the switching method, the ON/OFF ratio by which the electric power is supplied at a cycle which does not 
influence the motor speed is changed to perform a control. The switching method can be realized by periodically 
switching over Sy^ in Fig.26. The electronic circuit diagram close to this method is shown in Fig.27 and its functional 
model is shown in Fig.28. Figs.27 and 28 respectively show a model including the battery. 

[01 60] In Fig.27, a signal voltage E CNT of a rectangular wave applied from the outside is converted into a base current 
by a resistance R A of a transistor T R to make the transistor T R perform an ON/OFF operation. At this time, E 0 and R B 
respectively denote a battery voltage and an internal resistance. D,, in Fig.27 denotes a wheeling diode, and has the 
following two functions: One is a function of restoring the energy saved by the inductance of the motor winding generated 
when the transistor is turned OFF to the motor, and another is a function of guarding the transistor from the high voltage 
induced by the inductance. It is to be noted that a mechanism model of a switching signal generator is omitted. 
[0161] In the switching operation of Fig.28, a conductive state of ON is represented by a low resistance R ON and an 
interrupted state of OFF is represented by a high resistance Ro FF , in which both are switched over. Also, the function 
of the wheeling diode is represented by a condition judgment S m shown by a hexagon frame and a resistance R D of 
the conductive state. Namely, the condition judgment judges the voltage E CNT of the signal generator with the 
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detected voltage E^, the conductive resistance R ON is selected when 1 , and the interception resistance R OFF 
is selected when S W e=0\ Selection is performed by the switchover switch variables t and S WE _ 0 . 
[0162] Moreover, as for the accumulated energy of the inductance, the release period is judged by the counter elec- 
tromotive force of the inductance which has turned to be negative, so that the energy is restored to the motor through 
the conductive resistance R D of the wheeling diode. The functions of the braking resistance R L shown in Fig.26 and 
the resistance R D in Fig.28 appear to be the same at a glance. However, the former functions by a positive induced 
voltage of the motor, while the latter functions by a negative counter electromotive force by the inductance. The gov- 
ernment equation of the switching method taking these points into account is given as follows: 



0- 



^ON ^WE_1 + ^OFF ^WEO * 

1 ~ R D S WD 



.Eq.(44) 



[0163] It is to be noted that the switching signal generator is excluded. 
(3) Fixed current control method 

[0164] The electronic circuit for controlling the motor by continuously changing the motor current or voltage is as 
follows: 

[0165] The basic function of the motor is to convert the given current into the drive torque by the torque coefficient 
X T . Accordingly, when an electronic control is performed to the motor, it is a basic function of the fixed current control 
to freely set a current targeted and to maintain the current. Fig.29 shows one of control circuits for maintaining the 
output current l s in proportion to the value of a designated current l CNT . 

[0166] In Fig.29, a triangle shows an operational amplifier or the like which has a high amplification ratio inside. The 
current l CNT targeted is converted into the voltage by the resistance R s to be inputted to + side of the amplifier. Also, 
the controlled current l M is converted into the voltage by the current detecting resistance R E to be returned to - side of 
the amplifier, thereby performing a feedback control. It is to be noted that E 0 and R B are respectively the induced 
voltage and the internal resistance of the battery. 

[0167] Additional brief examination of the fixed current control will be described as follows: 
[0168] The voltages applied to the + and -sides of the amplifier are respectively V A , NP1 = RslcNT andV A_iNP2= R E , M» 
as seen from Fig.29. Also, assuming the amplification ratio of the amplifier is fi^p, the output voltage V A Jyjj is given 
as follows: 



V AOUT = ( R S I CNT _R E I m) A MP „ 

= V M+ R E I M ..JBq.(45) 

[0169] Also, the amplifier in Fig.29 has a high amplification ratio. Therefore, if equation (45) is arranged assuming 
1/A MP =0, the equation of the output current l M is given as follows, where R S /R E in the equation indicates the current 
amplification ratio: 

R s 

'm = rJi'cnt Eq.(46) 

[0170] V A 0UT of Eq.(45) cannot generate the voltage equal to or more than a supply voltage V B shown in Fig.29. 
Also, a voltage difference of an offset voltage V OFS is necessary between an output voltage upper limit value of the 
amplifier and a power supply voltage. This relationship is expressed by the following inequality: 

Vout*V v ofs Eq.(47) 

[01 71 ] The output voltage V A 0UT of the inequality (47) restricts the upper limit value of the output current l M , indicating 
the limit of the fixed current control. The output current l SL exceeding this limit is determined by the closed circuit of 
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the voltage of the battery induced voltage E 0 and the offset voltage E^g shown in Fig.29, as well as the voltage drop 
by the motor winding resistance R M , the current detecting resistance R E , and the battery internal resistance R B . How- 
ever, since E 0 is not included in the functional element of the fixed current control, it is necessary to replace E 0 with 
the supply voltage V B which can be added to the functional element. The current l SL considering this point is expressed 
5 by the following equation: 



10 



15 



Vd-E, 



B'^OFS 



Eq.(48) 



[0172] The current consumed by the internal circuit of the amplifier is very low. Therefore, by neglecting the current, 
equation (48) turns into the approximate equation of the current l B supplied by the battery. 

[0173] From the above-mentioned description, the fixed current control circuit in Fig.29 has a function of switching 
over the controlled current l M given in equation (46) to the uncontrolled current l SL given in Eq. (48) to be supplied. 
The functional model of the fixed current control in which the above-mentioned contents are arranged and modeled 
can be shown in Fig.30. 

[0174] Firstly, the government equation of the fixed current control can be obtained from Fig.30 as follows: 
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0 Swamp , Swamp o *cnt 
k e 

1 ^E^WAMP _J — S WAMp 



Eq.(49) 



[0175] In equation (49), S WAMP1 and S WAMP0 indicate switchover elements with the condition judgment forjudging 
the limit of the fixed current control. 
30 [0176] Finally, as for the fixed current control provides an in-controi while the current l B within the control range of 
Eq.(46) is equal to or less than the current l SL exceeding the control limit of Eq.(48), and a non-control when the current 
l B exceeds the current l SL . This condition judgment S WAMP provides the following equation. It is to be noted that 1 B and 
l SL in equation (50) are estimated state values. 
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if O b <i sl) then s wamp = 1 e,se s wamp = 0 



Eq.(50) 
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[0177] In Eq.(50), Sw^p = 1 means that the condition is fulfilled so that the current is directly supplied from the 
battery, while Sy^p^O means that the condition is not fulfilled so that the fixed current control is performed. 

3.2 Determination of initial characteristic value 
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[0178] Characteristics which have not yet been clarified are mainly transient characteristic values inside the motor. 
In order to understand the relationship between the motor function and the functional model, the characteristic value 
is initialized by providing the design value shown in the following Table 4 or an initial characteristic value obtained from 
a similar part, by applying a steplike voltage to the motor in the idle state, and by sequentially updating the characteristic 
value experimentally. It is a matter of course that the process and method performed here can be simulated by the 
computer 
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TABLE 4 



CHARACTERISTIC NAME 


SYMBOL 


UNIT 


CHARACTERISTIC VALUE 


INERTIA MOMENT 




[Nm/(rad/Sec 2 )] 


0.4x10*6 


INDUCTANCE 




[H] 


3.0x10-3 



[0179] Upon the identification, the simulation results at the time when the functional model, corresponding to model 
when T M = 0 in Fig.23, in the idle state is initialized by providing the characteristic shown in the Table 4 to Eq.(35) are 
shown in Figs.31 (1 ) and 31 (2). It is to be noted that for the steady internal characteristic value and the test condition 
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except those shown in the Table 4, the above-mentioned identification results of the steady state shown in the Table 
3 are used. Also, since the additional load T M is idle, it is made "0". 

[0180] Fig.31 (1 ) shows the motor current l M [A] and the angular velocity [rad/sec]. In Fig.31 (1 ), a thin dotted line 
indicates natural data of the measured current, a thick dotted line indicates a current l M MEG [A] where signal processing 
is performed to a measured result described later with a digital filter, a thick solid line"" indicates a current l M SIM [A] of 
the simulation result, and a thin solid line indicates the angular velocity [rad/sec] of the motor. 
[0181] Also, in Fig.31 (2), a thick solid line indicates a current deviation AI ERR [A] of an initial value obtained by 
subtracting the current l M MEG after being signal-processed with the filter from the current l M SIM determined by the 
simulation, and a thin dotted line indicates a current deviation AI P err [A] at the time when a moment of inertia J M 
described later is identified. 

[0182] In Fig.31 (2), the current l M shows a rise immediately after the power supply is turned ON where a high starting 
current flows, and decreases with a lapse of time, being stabilized to the idle current. Oppositely, the angular velocity 
u) M rises with a lapse of time to stabilize in the idle rotation number The change of the motor current indicates that the 
induced voltage of the motor is governed by the angular velocity through the motor constant M M . 

3.2.1 Noise process of measured current 

[0183] The motor current l M is measured by connecting a distributary resistance R s of 50 [mft] between the motor 
and the power supply. The starting current measured at this time is shown by a thin dotted line which varies extensively 
in Fig.31 (1). The extensive variance indicated by the actually measured value is due to the brush and a commutator 
of the motor, and is remarkably shown in a low rotation region immediately after the start. Particularly, the extensive 
variance which occurs at the rise of the rotation depends on the rotation angle of the commutator for the brush, which 
is remarkably shown in a small motor having a few contact pieces only of the commutator. 

[0184] The result of performing the signal processing to the measured data with the digital filter in order to remove 
the extensive current variance is shown by the thick dotted line in Fig.31 (1 ), for the filter characteristic, a low-pass filter 
with a cut-off frequency of 100 [Hz] is used. It is to be noted that since the phase delay of about 72 [msec] from the 
measured data has occurred due to filtering, the result of the signal processing in Fig.31 (1) is made coincident with 
the measured data by forwarding the rising point by the phase difference. Thus, the diagonally rising component having 
moved before 0 [sec] shown in Fig.31 (1 ) is neglected. 

[0185] The identification processing is performed by using the data to which the filter processing is performed. How- 
ever, it is expected that an error greatly influences the identification result since the current variance occurring at the 
start of the motor rotation is not completely excluded. Since this error is an uncertain variance dependent on a positional 
relationship between the brush and the commutator at the starting time described befo re, the identification is performed 
from this view point. 

3.2.2 Time history sensitivity of characteristic (process 38) 

[01 86] It has been clarified that all of the characteristic values which form the subject of the identification are included 
in the motor current l M obtained from the output state value of the output equation in the 5th row of Eq.(36). Therefore, 
in order to perform the identification, it is necessary to form a hypothesis by fully grasping the relationship between the 
characteristic values which form the subject and the motor current, and to verify it through the interpretation of the 
identification result. Accordingly, the followings are hypothesized for the influence of the characteristics on the motor. 

(j) It is supposed that the moment of inertia J M is proportional to the change of the angular velocity (angle accel- 
eration), provides strong influences during the acceleration when the change of the angular velocity is large, and 
does not provide influences in the steady state when the change is little. 

@ It is supposed that the inductance L M is proportional to the change of the motor current, provides influences 
when the current variance is large at the starting moment, and does not provide influences in the steady state 
when the current becomes stable. 

® It is supposed that the viscous resistance coefficient 5 M is proportional to the angular velocity, provides influ- 
ences as the load of the internal loss with the increase of the angular velocity, and assumes a fixed load in the 
steady state. 

@ It is supposed that the kinetic friction torque F MF is a constant and so is a characteristic not influenced by the 
state value change. The motor has a structure of reducing the friction torque as much as possible, so that the 
influence is little. 

[0187] It is seen from the above-mentioned contents that the respective ways in which the characteristic values 
contribute are different in the process of a time transition from the starting moment to the steady state. Accordingly, if 
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the sensitivity of the characteristic values for the change of the motor current l M is supposed to be the index of the 
identification, it can be easily presumed that the identification is possible. It is to be noted that the above-mentioned 
® and® are added for examining the time history sensitivity of the steady internal characteristic value. 
[01 88] The sensitivity of each characteristic value will be examined based on the above-mentioned result as follows 
(see flow chart in Fig.32): 

[0189] It has been described that the simulation result in Fig.31(1) has been performed at the initial values of the 
characteristic values, and the current deviation AI ERn at that time is shown in Fig.31 (2). In order to make this the index 
of the identification, the influence which the characteristic value variance provides to the current deviation AI ERR is 
clarified, so that the relationship between the characteristic value and the motor current can be grasped. This current 
deviation AI ERR can be expressed by the following equation (at step S3), with a current measured value l M MEG (at 
step S1 in Fig.33) and a simulation value l M SlM (at step S2) being made a deviation: 



[0190] If the deviation AI ERR lies within an allowable range, the transient identification is completed after verifying 
the identification result by the simulation described later (at steps S4 and S1 0). 

[0191] From equation (51 ), the value of the simulation is high with respect to the measured value when the current 
deviation Al ERR of Fig.31 (2) is positive, and otherwise low. Moreover, an extensive violent variance is included at the 
starting moment. Accordingly, in order to use the current deviation AI ERR as the sensitivity of the motor current l M for 
the characteristic change, it is required to remove the variance. Also, the offset quantity by the current deviation AI ERR 
is included in the sensitivity. 

[0192] A variance current deviation AI PERR of the characteristic representing the sensitivity is expressed by the 
following equation to remove the variance portion of the deviation: 



[0193] On the upper term of equation (52), the l P stM (at step S5) indicates the simulation value at the time when 
the transient internal characteristic value is changedTand the initial value current deviation AI ERR indicates the current 
deviation of the measured value I m _meg * or tne simulation value l M S | M of the initial value by Eq.(51). 
[0194] From this equation, the variance current deviation AI PERR of the characteristic value can use the deviation 
of the simulation value l P _ SIM of each characteristic value variance for the simulation value l M SIM of the initial value 
as a sensitivity as shown in the lower term of Eq.(52). However, since the identification accuracy is exceedingly eval- 
uated by the initial value current deviation AI ERR , it is necessary to proceed the identification while the AI ERR is contained 
within an appropriate range at any time during the identification. 

3.2.3 Identification method by time history sensitivity change 

[0195] Firstly, the result (at step S6) upon having determined the current deviation AI P ERR at the time when the 
characteristic values are changed by +1 0 [%] and -1 0[%] with the time history is shown in Figs. 33(2)-33(5). The four 
diagrams show the sensitivity seen with the time history at the time when the characteristic values are changed. Fig. 
33(1 ) is a diagram for comparing the sensitivities at the time indicated by thin dotted times in the lower four diagrams, 
so that the time history data are stored. 

[01 96] Also, in the respective current deviations of Fig.33(2)-33(5), J indicates the moment of inertia, L the inductance, 
D the viscous resistance coefficient, and F the kinetic friction torque, which are generally referred to "P (parameter)" 
as mentioned above. However, since the vertical axis scale in each diagram is the same as that in each diagram on 
the lower term, it is impossible to compare the four diagrams in parallel. Also, the solid line shows a result upon having 
increased the characteristic of the Table 4 by +10 [%], and the dotted line shows a result upon having decreased the 
same by -1 0 [%]. This calculation is performed to all of the characteristic values of the identification subject (at step S7). 
[0197] It is seen that the characteristic becomes linear, since the result upon having changed the characteristic value 
by ± 1 0 [%] is symmetrical with respect to positive and negative in Fig.33(1 ). The moment of inertia J M and the viscous 
resistance coefficient 6\, become directly proportional to the motor current, and the inductance L M becomes proportional 
to the motor current in which the positive and negative are inverted. These relationships almost agree to the above- 
mentioned hypothesis of the relationship between the motor current l M and the characteristic, so that the variance of 
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Eq.(51) 




...Eq.(52) 
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the deviation is removed. 

[01 98] The identification is performed with the current deviation AI P ERR at the time when the characteristic of Figs. 
33(1 )-33(5) is changed being made a time history sensitivity. In the identification, trials are repeated so that the current 
deviation lies within an allowable range by the following process (at step S4): 

[0199] Firstly, since the positive current deviation is large in the range of 0.02-0.06 [sec] from the initial value current 
deviation AI ERR of Fig.31(2), the moment of inertia J M in Fig.33(2) whose sensitivity is high in this region as shown in 
Fig.33(1 ) is selected (at step S8) to be reduced (at step S9) in the identification of the characteristic values. Since a 
variance value DI ERR at the starting moment comes near the upper side as shown in Fig.31 (2), the inductance L M of 
Fig.33(4) is increased. 

[0200] As for the order of identifying these characteristic values, the above-mentioned steps S2-S9 are sequentially 
executed from the value having larger coefficient with the mutual relationship between the initial value current deviation 
Al ERR and each variance current deviation AI PERR being made an index. Alternatively, when it is thought that the 
variances do not interfere with each other due to the time difference as shown in Figs.33{2) and 33(4), the identifications 
may be simultaneously performed by the steps S2-S9. 

[0201 ] It is to be noted that since the characteristic curve of the moment of inertia J M is similar to the current deviation 
Al ERR shown in Fig.31 (2) in the curve of the time history change, only the moment of inertia J M may be changed for 
the identification under the judgment that the current deviation AI ERR is most influenced by the moment of inertia J M . 
[0202] Among above-mentioned processes, the process firstly performed in which the updated value of the moment 
of inertia J M is determined is shown in Figs.34(1) and 34(2). 

[0203] Fig.34(1) is a diagram showing the current deviation AI ERR in Fig.31 (2) up to 0.06 [sec]. Point A in Fig.34(1) 
indicates a deviation quantity a which is desired to be reduced, and the position (point of 0.03 [sec]). Fig.34(2) shows 
the sensitivity of the moment of inertia J M at the point of 0.03 [sec] shown in Fig.33(1 ), in which the sensitivity point 0 
(abscissa) is adjusted to the point A in Fig.34(1). 

[0204] It is to be noted that the current deviation shown by the both ordinates is adjusted to the same scale. As for 
the deviation quantity a = 0.2 [A] which is desired to be removed from Fig.34(1), it is seen from Fig.34(2) that the 
sensitivity becomes negative. 

[0205] Accordingly, a correction ratio p = 16.8 [%] of moment of inertia J M is determined from an intersection B 
between the abscissa of the deviation 0 in Fig.34(1) and the extension of the negative sensitivity in Fig.34(2). From 
the Table 4, the updated moment of inertia J M assumes 0.4x 10' 6 x (1-0.168) = 0.3328 x 10' 6 . 
[0206] The result upon having updated the moment of inertia J M is a thin dotted line showing the variance current 
deviation Al P _ ERR of Fig.31 (2). The above-mentioned process is expressed by the following equations (at step S9): 



...Eq.(53) 



[0207] In equations (53), the upper term is an equation of the correction ratio p of the characteristic value, and the 
lower term is an equation of the characteristic value after the update. As for the variables in Eqs.(53), a(t) [A] is a 
current deviation quantity of the initial value current deviation AI ERR at a time T, Jt(t) [A/%] is the sensitivity of a char- 
acteristic value at a time T similarly, P M NEW is a characteristic value after the update, and P M OLD is a characteristic 
value before the update. 

[0208] When the motor current and the characteristic value have non-linear characteristics (time history sensitivity 
becomes asymmetrical) which are not proportional to each other, and the change range of the characteristic value 
greatly exceeds ± 10 [%], it is necessary to again determine the sensitivity from Eq.(52) to repeat the identification. 

3.3 Result of transient identification and simulation (at step S1 0) 

(1) Identification result and simulation thereof 

[0209] The characteristic values via the above-mentioned process of the transient identification are shown in Table 
5. As test data for identification, a result of an idle test is used in which the setting torque of the fixed torque transient 
test is assumed to be 0. 
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TABLE 5 



CHARACTERISTIC NAME 


SYMBOL 


UNIT 


CHARACTERISTIC VALUE 


INERTIA MOMENT 
INDUCTANCE 




[Nm/frad/sec 2 )] 
[H] 


0.335 x10" 6 
4.2x10-3 



[0210J The simulation result in which the transient identification result is reflected is shown in Fig.35. The simulation 
has been performed by providing the transient identification result of the Table 5 and the characteristic value of the 
steady identification result of the Table 3 to the government Eq.(35) for execution of the fixed torque transient test 
model shown in Fig.23. 

[0211] In Fig.35(1 ), the result out of consideration of the inductance L M is added by a thin dotted line and the natural 
data of the motor current are omitted, while the others are the same as Fig.31 . It is to be noted that the simulation out 
of consideration of the inductance L M has been performed by the same Eq.(35) with L M = 1 0* 12 [H] being made a minute 
value without influence. 

[0212] When both are compared, it is seen that by the influence of the inductance L M the starting current which 
steeply rises at the starting moment is saved as an energy to be released with the decrease of the current. 
[0213] As a result, the rising of the motor current and the peak value are suppressed in a solid line (L M > 0) with 
respect to the thin dotted line (L M = 0), and the falling part shifts to the right side. Also, the identification results are 
almost coincident as shown by the current deviation AI ERR of Fig.35(2), while the variance at the starting moment 
appears as it is in the form of the current deviation as mentioned before. 

(2) Simulation of actual load transient test model of motor 

[0214] Since this identified motor model becomes a virtual prototype which accurately reproduces a motor unit, a 
model of a virtual test which reproduces the driving state of the actual machine can be made by incorporating the 
functional model for operation in Fig.26 into the actual load transient test model shown in Fig.21 (see Fig.42 described 
later). The simulation result will be described: 

[0215] The simulation is performed by providing the characteristic value of the electrically driven arm of the following 
Table6toEq.(41). 



TABLE 6 



CHARACTERISTIC NAME 


SYMBOL 


UNIT 


CHARACTERISTIC VALUE 


DECELERATION RATIO 






1/200 


OPERATING ANGLE 




[rad] 


1.57 


STOPPER STIFFNESS 




[N/m] 


7.0x10* 


STOPPER RADIUS 


L S 


[m] 


0.04 



[0216] It is to be noted that the characteristic values other than those in the Table 6 are the same as those of the 
motor model after the identification. Also, the motor operation is performed by the power supply switch Sy^, the regular/ 
reversed rotation switch S^, and the braking switch S WB shown in Fig.26. 
[0217] Firstly, the simulation result of the motor start and stop is given as follows: 

[0218] The power supply switch S^ is turned ON to start the motor with the regular/reversed rotation switch S^ 
being made OFF and the braking switch S^ being made OFF. Then, the load of 2 [N • m] is provided after 120 [msec], 
and the power supply switch S ws is turned OFF after 280 [msec] with the load being provided. The simulation result 
is shown in Figs.36(1)-36(4). 

[0219] Figs.36(1)-36(4) sequentially show a motor current l M [A], a load torque T L [Nm], a motor voltage V M [V], a 
motor angular velocity [rad/sec], and an operation angle of arm e R [deg]. The motor current l M increases, and the 
angular velocity iom decreases from the moment the load is added. At this time, the change of the operation angle 9 R 
becomes rather slow. 

[0220] Also, when the power supply is OFF, the motor current l M assumes 0 [A], the angular velocity decreases, 
and the load drives the motor from about 31 0 [msec], so that the reversed rotation occurs. Also, as for the motor voltage 
V M after the power supply is turned OFF, the inductance L M generates a spike-like voltage at the moment the power 
supply is OFF, and then the induced voltage V w by the angular velocity of the motor occurs, so that the state where 
the motor is changed to a generator is described. 

[0221] The simulation result of the motor regular/reversed rotation, as well as the braking is as follows: 
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[0222] Figs.37(1 )-37(4) show an example in which the motor is stopped with and without the electric brake by OFF/ 
ON of the braking switch S WB . Figs.37(1)-37(4) show states as follows: The motor is started with the regular/reversed 
rotation switch S WL being made OFF and the power supply switch S ws being made ON to be compulsorily stopped 
with the stopper Then, the motor is rotated reversely in 0.6 [sec] after the regular/reversed switch S WL being made 
ON and is again compulsorily stopped with the stopper on the opposite side. After the lapse of 1.4 [sec], the power 
supply switch S^ s is made OFF, so that the motor is stopped. 

[0223] Figs.37(1 )-37(4) sequentially show the motor current l M [A], the operation angle 9 R [deg] of the arm, the motor 
angular velocity [rad/sec], and the motor torque T M [Nm]. Also, the solid line of Figs. 37(1 )-37(4) shows the state 
where the electric brake by the braking resistance of R L = 1 [ma] is applied, while the dotted line shows the state without 
the brake. Both of the regular and reversed rotations generate an attenuated oscillation at the time of collision at the 
stopper. 

[0224] Also, when the motor is compulsorily stopped by the stopper, a flexible structure of a stopper gum absorbs 
the maximum torque of the motor to be compressed. If the power supply is cut off in this state, the arm is pushed back 
by the compression power owing to the bend of the stopper. The effect of the electric brake appears at the stop state 
after 1 .4 [sec]. It is seen from the motor angular velocity and the operation angle 9 R of the arm that the motor idle 
by the stopper counterforce is well absorbed as shown by the solid line. 

4. Motor identification apparatus 

[0225] Hereinafter, an apparatus required for identifying the above-mentioned test model of the steady and transient 
states of the functional parts will be described. 

[0226] The functional parts incorporated into the products have a structure working with assembling the functional 
parts to the input/output portion and by mutually associating them. As shown in Fig.21 , the battery for driving, the 
electrically driven arm forming a load and the like are necessary e.g. for the motor test. Also, in order to control the 
motor, the operation system shown in Figs.25, 27, and 29 is necessary. 

[0227] The test model for thus performing an identification by using the test data of the actual machine requires 
modeling the functional parts which form the driving source and the load of the tested object. Subject to the influence 
by other functional parts in this way, it is possible to standardize the function and the behavior as the specification of 
the tested object and the test standard. If there is something which operates and reproduces the behavior just as 
standardized function and behavior, it can be used as a substitute common to each test method. 
[0228] If the functional model which forms the substitute of the driving source and the load is modeled and incorpo- 
rated into the test model as an identified model, the identification defined to the internal characteristic value of the 
tested object can be performed. Also, as for the testing unit and the test model, the group of the tested objects whose 
parts specifications are different can be replaced with block replacement means. 

[0229] Fig.38 shows an example in which the concepts are applied to the motor identification apparatus. In the 
identification apparatus of Fig.38, a motor model 1 has a block replacement structure for a test reproduction model 22. 
A reproduction model 2 is composed of a drive model 21 such as a battery and an operation model, a load model 22 
such as a load generator model, an observation model 23 for observing the state value of the motor model 1, and a 
characteristic update model 24 for updating the internal characteristic value. 

[0230] Also, a drive circuit unit 4 takes charge of the battery and the operation system which forms a driving source 
of a tested motor 3. A load generator 5 and an electric load 6 take charge of the mechanical load of the motor 3. The 
driving source and the mechanical load are modeled as a drive model 21 and the load model 22 already identified. 
[0231] Furthermore; as for the relationship between the actual machine testing unit and the models, the former drive 
model 21 has a structure of controlling the drive circuit unit 4 through a drive system control unit 71 , and the latter load 
model 22 has a structure of controlling the load generator 5 through a load system control unit 72. 
[0232] The measured value of the motor 3 is obtained at a measurement unit 8 to be provided to a calculation unit 
9. The calculation unit 9 receives the observation value of the observation model 23 and the measurement from the 
measurement unit 8 to update a characteristic update model 24. 

[0233] The motor drive model 21 shown in Fig.38 corresponds to the model in which the battery on the left side of 

Fig.21 and the operation system circuit of Figs.26, 28, and 30 are united. The load generator model 22 corresponds 

to the arm mechanism on the right side of Fig.21 and the fixed torque generation model of Fig.23. 

[0234] The operation and the behavior reproduced on the former model are reproduced at the drive circuit unit 4 by 

the drive system control unit 71 to be provided to the motor 3 of the tested object. Similarly, the load state provided to 

the motor 3 is reproduced on the load generator 5 by the load system control unit 72 controlling the electric load 6. 

[0235] Accordingly, if the motor functional model 1 shown in Fig.4(1) is connected to the reproduction model 2, the 

drive model 21 and the load model 22 preliminarily provided corresponding to the drive system and the load system 

turns into predetermined steady test models and control the drive circuit unit 4 and the electric load 6. 

[0236] Together with this operation, the steady test data generated by the motor 3 is measured by the measurement 
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unit 8, so that the calculation unit 9 which has received the data provides the government equation of the steady 
functional model i.e. the steady internal characteristic value to the motor model 1. Thus, the steady identification is 
completed. 

[0237) Then, the reproduction model 2 changes the drive model 21 and the load model 22 into the transient test 
model for the transient identification. Together with this operation, the transient test is executed, so that the transient 
test data are provided to the calculation unit 9 from the measurement unit 8. The calculation unit 9 compares the value 
obtained by substituting the already calculated steady internal characteristic value for the internal characteristic value 
inputted from the observation model 23 with the measured data from the measurement unit 8. The internal characteristic 
value of the motor model 1 is properly updated as mentioned above, by the characteristic update model 24 so that an 
error between them becomes small to lie within an allowable range. Thus, the transient identification ends, so that all 
of the identification works are completed. 

5. Virtual test 

[0238] The model which reproduces the function, the performance, and the characteristic values of the parts or the 
products composed of the parts is identified with the test data of the actuality, so that the internal characteristic value 
and the reproductivity are verified. Therefore, it can be utilized as a virtual prototype for performing various evaluation 
tests on the computer instead of the actual machine. An example where the virtual prototype is applied to the product 
development will be described. 

5.1 Relationship between product development and modeling 

[0239] The development process of a product and the formation process of a model are shown by the relationship 
in Fig.39. 

[0240] Fig.39 shows a relationship between a development process of a project stage, a design stage, a prototype 
and test stage, and modeling and identification performed in parallel with the development. Fig.39 shows a process 
including a concept of concurrent engineering in which development processes of the product are performed in parallel. 
[0241] On the other hand, the process of modeling has a flow as follows: An abstract product target at the project 
stage is firstly modeled, the model is sequentially detailed along with the development process, a validity of the project 
and design contents by the model is verified, and the identification is finally performed with the test data of the prototype, 
so that the virtual test is performed executing computer simulation with the identification model instead of the various 
evaluation tests. 

(1) Relationship between product development and modeling 

[0242] If this virtual test is applied to the development of a car, for example, the process is given as follows: 
[0243] As for the abstract model at the project stage such as a model reproducing a car performance which forms 
the target characteristics of the projected car such as an acceleration time of 0-400 [m], a stationary fuel cost perform- 
ance, and an acceleration performance, a possibility of achieving the target characteristics is verified by the model 
assembled with the characteristics such as weight, resistance, power characteristic, and gear ratio. 
[0244] In the next product part formation and basic design, the model is disassembled following the formation method 
of the basic functional parts such as an engine, a transmission, a body realizing the target characteristics, and modeling 
is performed with the internal characteristic values of the functional parts being clarified, so that the function and the 
performance required for the functional parts, and the possibility of realizing the target characteristic value are verified. 
[0245] Similarly, in the detailed design, the modeling is detailed according to the contents of the detailed design of 
an engine, a clutch, a brake, a transmission, a control unit, and the like, whereby the validity of the design contents in 
the details such as torque variance, transmission shock, vibration noise, and the control algorithm is verified with the 
model. 

[0246] The model made by the above-mentioned process is sequentially disassembled into the lower model from 
the abstract model at the project stage along the part formation, so that a systematic model hierarchized from the 
product to the part and its details is formed. 

[0247] Also, in the above-mentioned process of the project, design and modeling, the development is proceeded by 
the design drawing of the actuality realizing the target characteristic, and the model verifying the design drawing. 
However, both are eventually hypotheses. Accordingly, it is necessary to test both by an experimental car made based 
on the design drawing or the like, to identify the model, and to verify the validity of the hypotheses. 
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(2) Relationship between product development and identification 

[0248] By the above-mentioned reason, the identification shown in the lower side of Fig.39 is performed based on 
the data of the functional parts and a car test. The identification and the test necessary therefor are performed separately 
in the above-mentioned steady and transient states. The model is converted in accordance with the test method as 
shown in Figs.9-1 3 for the steady test, and in Rgs.21 -30 for the transient test, so that the identification of the test model 
is performed by the test data. In the identification performed here, the identification result of the test models firstly 
performed per functional part are integrated along the part formation to be made the identification of the product model. 
[0249] It is supposed that the functional parts of the battery, the motor, and the electrically driven arm in Fig.21 are 
individually modeled, for example, and the product in which the parts are mutually connected is an electric actuator. 
The model of the electric actuator can be represented by a hierarchized model where the models of the battery, the 
motor, and the electrically driven arm are made lower level. Also, it is seen that the equation in which the models are 
integrated turns into Eq.(34). 

[0250J Accordingly, if the motor at the center of Fig.21 and the electrically driven arm on the right side of Fig.21 are 
individually identified, and they are made the characteristic value of Fig.21 , the electric actuator on the upper hierarchy 
is naturally identified. The model which can faithfully reproduce the characteristic value and the performance of the 
product along the mode of the actuality is called a "virtual prototype". 

[0251 J As an applying method of this hierarchized model, it is possible to convert the electric actuator model into the 
test model when the battery and the electrically driven arm are already identified in Fig.21 , to make the identified internal 
characteristic values fixed values, and to perform a partial identification defined to an unfixed internal characteristic 
value of the motor model. 

5.2 virtual test by virtual prototype 

[0252J Thus identified part model is confirmed as a virtual prototype of the part reproducing the function and the 
performance held by each part. Similarly, the product model in which the part models already identified are mutually 
connected and integrated means that the model is confirmed as a virtual prototype reproducing the product. Also, these 
models result in the identification by converting the models into test models in accordance with the above-mentioned 
test methods, comparing and evaluating the reproducibility and the test result of the actuality, and verifying the difference 
between both. 

[0253] Accordingly, the virtual prototype of the identified product and parts can be evaluated as to the validity by 
performing the computer simulation instead of the various tests of experimental articles and the like. Performing the 
evaluation test of the actuality on the computer by utilizing the virtual prototype is called a "virtual test" of the parts or 
the product. The simulation which has verified the validity of the identification result obtained in Figs.37-38 is also one 
example. 

(1 ) Internal formation of virtual test model 

[0254] The virtual test is usually applied with the test standard determining the same environment condition as the 
evaluation test of the actuality, a test condition and performance such as driving operation, and an evaluation standard 
such as a success/failure judgment, a performance prediction, a presence/absence of abnormal phenomena, and the 
like by the reproduction data of the performance and the behavior. This can be typically shown in Fig.40. 
[0255] Fig.40 is a basic arrangement of a virtual test model 1 0 where a virtual prototype 1 1 and a test standard group 
12 applied thereto are modeled to be incorporated. The virtual prototype in this arrangement comprises two kinds of 
functional part models, an already identified functional part model and a functional part model in which design values 
and characteristic values of similar parts are provided upon the virtual test. 

[0256] Both are distinctive in that the former is used as a model of a part having a basic function governing the 
function and the performance of the product, while the latter is used as a model extended for various products whose 
characteristic values are different in the application, the structure, and the shape of the products. The product model 
having both models united indicates that the evaluation test of various products and the test evaluation of various 
product groups where the derived functional models (latter) of different characteristic values are combined with the 
identified basic functional model (former) can be performed by the virtual test. 

(2) Virtual prototype example of virtual testing car 

[0257] For example, a simple model of the virtual prototype for a car in Fig.40 can be shown in Fig.41 . Fig.41 is an 
example of a virtual prototype showing a functional model where the functional parts of the engine and the transmission 
are modeled to be integrated as a power train (P/T), and a functional model where various types of cars which are 
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different in their structure and the like are modeled as a body, and both are integrated as a car. Between the models 
of the functional parts, the engine and the transmission are related with the angular velocity oo e that is a potential 
quantity, and with the torque T e that is a flow quantity, while the transmission and the body are related with the torque 
v b that is a potential quantity, and with the torque f b that is a flow quantity. 
5 [0258] The characteristics in the models are as follows: 

[0259] The characteristics of the engine are a moment of inertia J e> a viscous resistance coefficient 8 e , and an internal 
torque T e - 8 e and T e among them respectively indicate a partially linearized gradient, and a generation torque when 
the engine rotation number is "0", in the rotation/torque characteristic in which a throttle divergence ot e , an atmospheric 
pressure P e , and an air temperature tg are inputted. 
w [0260] The characteristics of the transmission are the gear ratio N m , the final gear ratio N d , and the stiffness of an 
output axis C p . Finally, the characteristics of the body are the tire radius R T the mass of the body M b , the resistance 
during running D b , and the ascent of a running road 6. Since the characteristics of the body, R T M b , and D b are different 
according to the body structure or the like, they are required to be replaced according to the type of car. 
[0261 ] Moreover, the control characteristic values of the engine control unit and the transmission control unit respec- 
ts tively controlling the engine and the transmission are preliminarily identified by the test data. 

[0262] In order to perform a virtual test under the same condition as the car driving state by using the virtual prototype 
in Fig.41 , the following model is further required. Firstly, for the gradient 6, the atmospheric pressure p e , the air tem- 
perature t e and the like, the environment condition model provided as an environment condition for using cars is re- 
quired. Furthermore, for the throttle divergence cc e , and the gear ratio N g of the transmission, the operation model which 
20 provides the driving operation condition of cars is required. 

[0263] Finally, in order to evaluate the reproduction result of the virtual prototype in which the virtual test has been 
performed, the observation model for observing a> e , V b , T e , f b , and the like is required. The model for the conditions 
and for performing the evaluation is modeled based on the evaluation test methods of the test standard models shown 
in Fig.40. 

25 [0264] The virtual test performed by the model can be applied to various types of cars by incorporating the charac- 
teristic values of the body model which are different depending on the type of car into the virtual prototype of the 
functional part in the P/T which faithfully reproduces the characteristic value, performance, and behavior of the actual 
machine. The virtual test using the virtual prototype of the car model can be substituted for the various tests performed 
to various types of car on which the same P/T is mounted. 

30 [0265] In the development of cars so far, the evaluation tests extending over a wide range have been performed to 
the cars having many body structures, so that the validity as a commodity has been verified. In such vast numbers of 
tests, the number of cars made on the experimental basis can be reduced and the development period required for 
the test can be shortened by the virtual test performed with the body characteristic value in the virtual prototype being 
replaced and with modeled test standard groups being combined. As a result, cost for a development can be reduced. 

35 

5.3 Virtual testing system (model) 

[0266] The basic arrangement of the actual machine test and the identification as well as the virtual test by the virtual 
prototype can be shown as in Fig. 42. 
40 [0267] It has been already described for the virtual test of the virtual prototype that the tested model which forms a 
subject of the actual machine test can be converted into the test model, the test standard can be modeled as the driving 
operation model and the environment condition model, and the execution result of the virtual prototype can be modeled 
as the observation model. 

[0268] It is possible to compare and evaluate both of the observation data reproduced by the virtual test and the test 
45 data of the actual machine test and to perform the identification. This indicates that the virtual test model and the test 
model for the identification can be made common. 

[0269] Namely, the virtual testing system (model) can have the consistency of the actual machine test, the identifi- 
cation, and the virtual test by incorporating an identification model 9 (corresponding to a calculation unit in Fig.38) 
which models the identification technique and process to update the internal characteristic value of the virtual prototype, 

so and an evaluation model 35 in which the virtual test verifies the validity of the reproduction data. 

[0270] In Fig.42, the actual machine test data are inputted to the evaluation model and to the identification model, 
and the result of the evaluation model is inputted to be related. For both models, the following processes are performed: 
In case a problem arises in the result of the success/failure judgment in the virtual test, the former performs the actual 
machine test to evaluate again as shown in a dotted line, and in case there is a problem in the identification result of 

ss the virtual prototype, the latter can perform the identification again. 

[0271] Main contents of each portion in the virtual testing system shown in Fig.42 are given as follows: 

© The virtual prototype 1 corresponds to the motor model in Fig.38, to a model which reproduces the function, 
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the performance, and the characteristic value of the product and part, which is composed of a mechanical system 
model relating to the mechanism and the structure, and to the control system model for controlling them. 
(D A driving operation model 31 is a model for operating the virtual prototype 1 by the same driving condition as 
the actual machine. For example, there are an accelerator operation, a transmission operation, a braking operation 
5 of a car, and the like. 

(3) An environment condition model 32 is a model for reproducing the influence of the environment with the virtual 
prototype 1 being made operate under the same condition as the environment for using the actual machine. For 
example, there are a running road, a temperature, a humidity, an atmospheric pressure, and the like when the car 
makes a market running. 

10 @ An observation model 23 is a model for reproducing the data, on the model, at each measurement point meas- 
ured in the actual machine test as shown in Fig.38, so that the data which are difficult to be measured in an actual 
machine can be observed on the model. 

© An internal characteristic update model 24 is a model for updating the internal characteristic value of the virtual 
prototype based on the identification result, as shown in Fig.38. 
*5 (D An actual machine test portion 33 performs a test which has been performed in the prior art, and a measurement 
unit 34 in the actual machine test portion 33 shows a unit for taking the test data into a computer on real time basis. 

[0272] The above-mentioned operation can be summed up as follows: If the driving operation condition and the 
environment condition are provided to the virtual prototype 1 from the models 31 and 32 respectively, the virtual pro- 
20 totype 1 provides the reproduction data obtained as a result of the simulation at this time to the evaluation model 35 
from the observation model 23. The evaluation model 35 evaluates, instead of the actual machine test, whether or not 
the reproduction data of the simulation result is valid (validity of product/part) by comparing the evaluation reference 
determined by the test standard with the reproduction data. 

[0273] In case of a failed evaluation, the evaluation model 35 instructs the virtual test portion 30 to change the driving 
25 operation condition and the environment condition, and to conduct the identification again if the deviation between the 
actual machine test data and the reproduction data is large. Thus, as shown in the dotted line in Fig.42, the driving 
operation condition and the environment condition are respectively provided to the actual machine test portion 33 from 
the models 31 and 32. 

[0274] Accordingly, if the deviation between the actual machine test data from the actual machine test portion 33 
30 through the measurement unit 34 and the reproduction data is large, the evaluation model 35 instructs the identification 
model 9 to conduct the identification again. Since the identification model 9 corrects the internal characteristic value 
of the virtual prototype 1 by the internal characteristic update model 34 after receiving the actual machine test data, 
the evaluation is performed again at the evaluation model 35 through the observation model 23. As a result, the iden- 
tification model 9 updates the internal characteristic value of the virtual prototype 1 so that the deviation between the 
35 actual machine test data and the reproduction data becomes minimum. 

[0275] When the evaluation ends, another modification (e.g. modification according to a type of car) in the virtual 
prototype 1 is evaluated. Particularly, by performing the identification to the internal characteristic value (control pa- 
rameter) of the control model in the virtual prototype, the virtual prototype 1 can be tuned, according to the type of car, 
to an optimum state. 

40 

5.4 Actual machine test of transmission applying virtual test 

[0276] In case the actual machine test is performed to a unit of the functional part of the tested object, the functional 
part of the drive system for driving the part and the functional part of the load system which forms the load are required. 
« The method in which the virtual prototype, where such a drive system functional part and a load system functional part 
are identified, is incorporated into the testing unit and the actual machine test of the tested object is performed, will 
now be examined. 

[0277] Fig.43 shows an example of a testing unit of a transmission included in the virtual prototype. 

[0278] Fig.43 shows a unit for reproducing the drive characteristic of the engine and the load characteristic of the 

50 car by the drive motor and the load generator, and for performing the actual machine test to the transmission proper 
along the actual use state. This unit derives the performance required for the running of the car from the engine per- 
formance which forms the drive source, and the evaluate test about various performances relating to the acceleration 
such as an acceleration performance or a gearshift shock is performed by the actual machine. 
[0279] In Fig.43, the testing unit 33 is composed of the drive motor 4 for driving the transmission, the load generator 

55 5 which forms the load, the electric load 6, and the measurement unit 8 for measuring the driving state. Also, to the 
testing unit 33 the virtual prototype 1 1 of the drive motor 4, the load generator 5 and the transmission 50 is connected 
for reproducing the driving state of the car. The testing unit 33 comprises the drive system control unit 71 and the load 
system control unit 72 which faithfully reproduce the result at the drive motor 4 and the load generator 5. 
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[0280] These control units 71 and 72 provide the drive torque of the engine to the transmission 50 by the drive motor 
4, and reproduce the load torque provided to the transmission 50 through the body by the load generator 5 and the 
electric load 6. Both models of the drive motor 4 and the load generator 5 can reproduce the function, the performance, 
and the behavior of the engine and the body in the testing unit 33 by the identified virtual prototype 11 . 
[0281] The virtual prototype in which these drive motor 4, transmission 50 and load generator 5 are modeled is shown 
in Fig.44. It is to be noted that the operation system models 31 , 38, 39, and the environment system model 32 inputted 
to the transmission model 1 in Fig.43 are the models of the test standard mentioned above. The former represents the 
model which inputs the driving operation condition of the driver, and the latter inputs the environment condition such 
as an air temperature which influences the internal characteristic value of the transmission. Also, the observation 
system models 23, 37, and 40 are models for observing the state values inside the virtual prototypes. 
[0282] The motor model 1 3 of the drive system in Fig.44 is composed of the models of the motor, the battery, and 
the fixed current control, respectively shown in Figs.6, 4A, and 30. It is to be noted that a part of the internal characteristic 
values of the motor model 13 are omitted for simplified description of the model. Similarly, on the load side, the load 
generator model 1 4 is composed of the braking resistance in Fig.4(1 ) and the electric load is composed of the braking 
resistance on the left side of Fig.13, where a part of the internal characteristic values is omitted. 
[0283] In Fig.44, the angular velocity <o e and the input torque T e connecting the drive system motor model 13 and 
the transmission model 1 are provided as target values of the drive system control unit 71 . Similarly, a velocity V b and 
a drive force f b of the load generator model 14 are provided as target values of the load system control unit 72. Also, 
the drive torque T e of the drive motor 4 can be reproduced by controlling the control current ! CNT with the operation 
system from the fixed current control model in Fig.30 and by adjusting the current of the drive motor. 
[0284] Also, the load torque of the load generator 5 can be reproduced by being changed into the braking resistance 
value through the operation system from the motor model of the braking test in Fig. 1 3. Also, as for the load torque, the 
fixed current control model and its unit in the same Fig.30 as the drive motor may be incorporated instead of the braking 
resistance. 

[0285] Moreover, into the testing unit 33 shown in Fig.43, the internal characteristic update model 24 is incorporated 
in which the evaluation model (see Fig.42) provided in the calculation unit 9 evaluate a difference if generated between 
the test result of the transmission 50 and the reproduction result of the virtual prototype 11 , and the identification model 
(see Fig.42) similarly provided in the calculation unit 9 identifies the internal characteristic value in the transmission 
model (virtual prototype) 1 . 

5.5 Actual machine test example of power train (P/T) applying virtual test 

[0286] An example will be examined in which the actual machine test of the P/T where the engine and the transmis- 
sion, except the motor 4 on the drive side and the virtual prototype 1 3 and the drive system control unit 71 leading to 
the motor 4, are united is performed in the actual testing unit shown in Fig.43. 

[0287] The example of the actual machine testing unit having the P/T and the load generator is shown in Fig.45. 
[0288] In this example, the same P/T 51 is mounted on the bodies whose structures are different shown in Fig.41 , 
and the actual machine test is performed to the various cars, so that the adaptability of the P/T 51 for the cars having 
the different body characteristic values is verified. 

[0289] The virtual prototype 1 of the P/T 51 requires many kinds in which transmissions such as an automatic trans- 
mission, a manual transmission, and a continuously variable gear are combined in addition to kinds of motors such as 
a motor, a gasoline engine, and a diesel engine. Accordingly, these virtual prototypes 1 are replaced with a model 
replacement unit 25. 

[0290] Also, the virtual prototype 1 of the P/T 51 has a control model for reproducing the software of the control unit 
of the engine and the transmission, so that the virtual prototype 1 can control the model of the virtual prototype 1 and 
can directly control the actual machine of the engine and the transmission through a control input/output unit 36. 
[0291] Accordingly, when adaptability of the body performance reproduced at the actual machine P/T 51 and the 
load generator 5 is worsened, the internal characteristic value of the control model can be changed through an internal 
characteristic update unit 24. Then, the measure effect is firstly confirmed in the virtual prototype 1 , and secondly the 
adaptability is again evaluated by the actual machine test to be confirmed. 

[0292] In case a low-speed region torque of engine mounted on a car whose body weight or running resistance is 
large is insufficient and the starting acceleration power is insufficient, the control characteristic value of the control 
model is updated through the internal characteristic update unit 24 in order to improve the acceleration performance 
and increase the fuel amount at starting the acceleration. The measure effect of the control amount change and the 
bad influence on the fuel-efficiency and the exhausted gas worried because of the increased amount are firstly con- 
firmed by the virtual test in the simulation of the virtual prototype, and the actual machine evaluation test of the P/T 51 
is finally performed by using the virtual testing unit. 

[0293] By this process, the tests for reproducing bad influence which examines the measure effect and many test 
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items can be verified by the simulation of the virtual prototype in a short time, so that it becomes possible to define the 
actual test items requiring a long time based on the verification result, and to perform an efficient test operation. 
[0294] As described above, a characteristic value identification method and an apparatus therefor according to the 
present invention is arranged so that a functional model of a part is prepared based on a potential quantity and a flow 
quantity representing a strength and a quantity of energy applied to the part, a steady internal characteristic value of 
the functional model in a steady state is identified, and a transient internal characteristic value of the functional model 
in a transient state is identified by using the identified steady internal characteristic value. Therefore, the function of 
all of the articles which govern two dimensions of the potential quantity and the flow quantity prescribing the energy 
can be modeled. 

[0295] Also, since the transient identification is performed to the internal characteristic value of the model after the 
steady identification, the steady identification without an interference of a transient state can be performed. Since the 
transient identification is performed based on thesteady identification, there is an effect of obtaining an accurate internal 
characteristic value by a simplified identification process. Also, it is possible to faithfully reproduce products/parts in 
the steady state and the transient state. 

[0296] Moreover, by manufacturing an apparatus which uses such a characteristic value identification method, it 
becomes possible to promptly perform the identification of the functional model of the same kind. 
[0297] Furthermore, the characteristic value identification apparatus according to the present invention is arranged 
such that the functional model having the characteristic value identified by the identification apparatus is incorporated 
in a virtual testing system as a virtual prototype, reproduction data are obtained from the virtual prototype by providing 
a driving operation condition and an environment condition, actual machine test data obtained by the driving operation 
condition and the environment condition are compared with the reproduction data, and a re-identification is performed 
depending on the comparison result as the occasion demands. Therefore, it becomes possible to omit the actual 
machine test performed in a development process of a design, prototyping, and a test, to shorten a period, and to 
reduce cost for a deveiopment. 



Claims 

1 . A characteristic value identification method comprising: 

a first process for preparing a functional model of a part based on a potential quantity and a flow quantity 
representing energy applied to the part, 

a second process for converting the functional model into a steady functional model in a steady state to identify 
a steady internal characteristic value, and 

a third process for identifying a transient internal characteristic value of the functional model in a transient 
state by using the steady internal characteristic value. 

2. The characteristic value identification method as claimed in claim 1 wherein the second process includes; 

a first step for determining an internal characteristic value of at least one steady test model from the steady 
functional model, 

a second step for collecting steady test data by performing a test corresponding to the steady test model, and 
a third step for identifying a steady internal characteristic value of the internal characteristic value based on 
the steady test data. 

3. The characteristic value identification method as claimed in claim 2 wherein the first step determines the internal 
characteristic value from a government equation in the steady state of the functional model. 

4. The characteristic value identification method as claimed in claim 3 wherein the third step converts the government 
equation into a recurrence equation to determine the steady internal characteristic value from a recurrence coef- 
ficient of the recurrence equation. 

5. The characteristic value identification meth od as claimed in claim 2 wherein the third step divides the steady internal 
characteristic value into a known factor and an unknown factor to identify the steady internal characteristic value 
of the unknown factor. 

6. The characteristic value identification method as claimed in any one of claims 1 to 5 wherein the third process 
includes: 
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a first step for determining an internal characteristic value of at least one transient test model in a transient 
state of the functional model, 

a second step for collecting transient test data by performing a test corresponding to the transient test model, 
a third step for applying the steady internal characteristic value to the internal characteristic value of the tran- 
sient test model to generate transient phenomenon reproduction data, and 

a fourth step for correcting the transient phenomenon reproduction data based on an error between the tran- 
sient phenomenon reproduction data and the transient test data, thereby identifying a transient internal char- 
acteristic value. 

7. The characteristic value identification method as claimed in claim 6 wherein when the error does not lie within an 
allowable range the fourth step repeatedly corrects a predetermined transient internal characteristic value within 
the transient phenomenon reproduction data until the error lies within the allowable range, and determines the 
transient internal characteristic value to be identified when the error lies within the ailowable range. 

8. The characteristic value identification method as claimed in claim 7 wherein the fourth step preliminarily calculates 
a variance deviation, as a time history sensitivity, to an initial value at a time when each transient internal charac- 
teristic value is increased or decreased at a fixed ratio, and selects a transient internal characteristic value having 
a maximum sensitivity within the time history sensitivity as the predetermined transient internal characteristic value. 

9. The characteristic value identification method as claimed in claim 7 wherein the fourth step preliminarily calculates 
a variance deviation, as a time history sensitivity, to an initial value at a time when each transient internal charac- 
teristic value is increased or decreased at a fixed ratio, and selects a transient internal characteristic value having 
the time history sensitivity similar to the error as the predetermined transient internal characteristic value. 

10. The characteristic value identification method as claimed in claim 9 wherein the fourth step simultaneously selects 
a plurality of transient internal characteristic values having different maximum sensitivity times as the predeter- 
mined transient internal characteristic value. 

11. A characteristic value identification apparatus comprising: 

block replacement means for a functional model of a part prepared by a potential quantity and a flow quantity 
representing a strength and a quantity of energy applied to the part, 

test reproduction means for reproducing at least one steady test model in a steady state of the functional 
model and at least one transient test mode! in a transient state, 

testing means of the part for performing a steady test and a transient test respectively corresponding to the 
steady test model and the transient test model, 

measurement means for collecting steady test data and transient test data at a time when a steady test and 
a transient test of the part are performed by the testing means, and 

calculating means for identifying a steady internal characteristic value of the steady test model by using the 
steady test data, for applying the steady internal characteristic value to the transient test model to generate 
transient phenomenon reproduction data, and for correcting the transient phenomenon reproduction data 
based on an error between the transient phenomenon reproduction data and the transient test data, thereby 
identifying a transient internal characteristic value. 

12. The characteristic value identification apparatus as claimed in claim 11 wherein when the error does not lie within 
an allowable range the calculating means repeatedly correct a predetermined transient internal characteristic value 
within the transient phenomenon reproduction data until the error lies within the allowable range, and determine 
the transient internal characteristic value to be identified when the error lies within the allowable range. 

13. The characteristic value identification apparatus as claimed in claim 11 wherein the calculating means preliminarily 
calculate a variance deviation, as a time history sensitivity, to an initial value at a time when each transient internal 
characteristic value is increased or decreased at a fixed ratio, and select a transient internal characteristic value 
having a maximum sensitivity within the time history sensitivity as the predetermined transient internal character- 
istic value. 

14. The characteristic value identification apparatus as claimed in claim 11 wherein the calculating means preliminarily 
calculate a variance deviation, as a time history sensitivity, to an initial value at a time when each transient internal 
characteristic value is increased or decreased at a fixed ratio, and select a transient internal characteristic value 
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having the time history sensitivity similar to the error as the predetermined transient internal characteristic value. 

15. The characteristic value identification apparatus as claimed in claim 1 3 wherein the calculating means simultane- 
ously select a plurality of transient internal characteristic values having a different maximum sensitivity time as the 
predetermined transient internal characteristic value. 

16. A virtual testing system which incorporates a functional model, as a virtual prototype, having an internal charac- 
teristic value identified by a characteristic value identification apparatus claimed in claim 11 comprising: 

condition assigning means for assigning a driving operation condition and an environment condition to the 
characteristic value identification apparatus, 

observation means for observing reproduction data obtained by the virtual prototype when the driving operation 

condition and the environment condition are assigned, and 

evaluation means for evaluating an observation result of the observation means. 

17. The virtual testing system as claimed in claim 16, further comprising another measurement means for measuring 
actual machine test data at a time when the driving operation condition and the environment condition are provided 
to an actual machine which forms a subject of the virtual prototype, and 

re-identification means of the virtual prototype, 

the evaluation means comparing an output of the measurement means and the observation result, and mak- 
ing the re-identification means re-identify the virtual prototype according to the comparison result. 

18. The virtual testing system as claimed in claim 1 7 wherein a fixed virtual prototype is incorporated into a part of a 
drive system and a load system connected to the part as the virtual prototype, the testing means perform a test 
corresponding to the fixed virtual prototype, and the evaluation means at this time make the re-identification means 
perform a re-identification according to the comparison result. 
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FIG.2 
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FIG.4 
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FIG.15 
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